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XXVI. Theory of Elasticity in General Coordinates. 


By A. E. GREEN and W. Zerwna, 
Department of Mathematics, King’s College, Newcastle-upon-Tyne f. ° 


[Received December 22, 1949.] 


$1. Inrropuction. 


‘THE tensor calculus (Ricci-calculus) has for many years been one of the 
standard tools in certain branches of mathematics and physics as, for 
example, in the general theory of relativity. So far, however, only a few 
writers appear to have used. tensor calculus for the presentation of a 
general theory of elasticity, either for infinitesimal or finite displacements. 
Elasticity theory is often developed in special rectangular cartesian 
coordinates and, particularly for finite displacements, the algebraic 
difficulties are considerable. Also, transformation to any other coordinate 
system is usually difficult and laborious. The tensor treatment avoids 
these difficulties and is in many ways easier to follow once the elementary 
ideas of tensor calculus have been mastered. 

In 1937 an important paper was published by Murnaghan (1937), in 
which he obtained some of the basic relations of the theory of elasticity 
for finite displacements J, using general tensor notations. This seems to 
be a fundamental paper in many respects but the present writers believe 
that the theory developed there is not yet in its simplest form. Since 
1937 papers have been published § in which tensor calculus notations of 
various types are used for elasticity theory, but the notations and the 
final forms of some of the results differ considerably, and are not always 
very simple. 

An examination of previous work suggests that this is an appropriate 
moment in which to attempt a simplification of the presentation of the 
theory of elasticity in general coordinates. While acknowledging the help 
derived from previous papers the present writers seek to give a reasonably 
complete account of elasticity theory in tensor notation, for finite dis- 
placements, in which it is believed that the treatment is straightforward, 
and that the methods used and the relations which are obtained contain 
many new features. 


g Communicated by the Authors. 
+ These relations have recently been expressed in vector and dyadic notations 


py “Milne- Thomson (1948). 
§ Some of these papers are given in the references. 
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§ 2. DispLACEMENT VECTOR. DEFORMATION. 


Let every point of a continuous three dimensional body, called 
briefly the body Bg, be at rest, at time t=0, relative to a fixed rectangular 
Cartesian system of axes a; whose origin is 0. The position vector of a 
typical point P, of the body B,, referred to o as origin, has the form 


r=f(7,. Vo, ali . . . . . . . . (2.1) 
and the line element is given by 
dse=dax' dx' ee ad itr O48 


\ - 
where the usual summation convention is used and the latin indices take 
the values 1, 2, 3. 


0 


We suppose that the positions of the points of the body B, change 
so that at time ¢ a typical point Py has moved to P. The point P may be 
referred to a new fixed rectangular Cartesian system of axes y; whose origin 
is at.O, and the position vector of P referred to O is (see fig. 1) 


R=R(¥1, Y2, Ys)- sO ie eke Shi cS 


The position vector v of the point P relative to P) represents the 
displacement of Py when the points of By move from hae positions ab 


time t=0 to their positions at time f¢. We therefore call 
v=R-r-+c, eer! my (2.4) 


the displacement vector, where the constant vector ¢ defines the position. 
of the origin O relative to the origin o. 
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Except for the assumption that the body B, is continuous no other 
assumptions have been made about the material of the body or about the 
way in which the motion takes place. We now make some assumptions 
about the way each point P, of By moves to its corresponding position 
P at time ¢. This will imply that the material of the body has certain 
physical properties. We shall not, however, discuss the nature of these 
physical properties. Although we are guided to some extent by physical 
ideas in making our assumptions, the justification for applying the theory 
_ to physical bodies must rest on the comparison of results obtained from 
theory with those found by experiment. 

We wish to study motions of points of the body By which take place 
continuously and which are such that any curve in By moves so that it 
always consists of the same particles. In addition, any line element ds, 
of By always remains of the same order of magnitude. We therefore 
assume that each point P, at time ¢, is related to its original position P,,, 
at time t=0, by the equations 


Y = Y (X1, Lg, Xz, t), SALE sa ana, Bon P45 | 


where y; are continuous single-valued functions of #,, 7, 7; and ¢, and 
where, for each ¢, these relations are reversible, so that 


Li=L(Y 15 Yos Yao 4), Series eh ea et. i (2.6) 


a, being single-valued functions of y,, ¥2, y3 for each value of t. 
We also assume that, at a given time ¢, the differentials dy’ exist and 
are related to the differentials dx’ by the equations 


= OYE ayk 0% ak 
de ep line : Sr amg Se pe Mies vee” (aie ) 
These assumptions impose certain restrictions on the partial derivatives 
of y; with respect to #,. Additional restrictions will also be necessary in 
later parts of the theory but we shall generally suppose, without further 
discussion, that all conditions which are necessary for the mathematical 
validity of the theory are satisfied. 

If, as a result of the motion of the points of Bo, the relative positions of 
neighbouring points of the body are altered, so that the length dsq of the 
line element changes to ds, the body B, is said to be strained or deformed.. 
The strained or deformed state of the body, at time t, is denoted by B.. 
The square of the line element in the strained body B may then be written 


ds?=dy' dy’, CI a ee eg RS) 


where the differentials dy’ are to be evaluated for a fixed time ¢ and are 
related to dx’ by equation (2.7). If ds=ds, for all points and for every 
time ¢ the body is said to be rigid. 
The original body B, may be described by a general curvilinear set of 
coordinates 0, so that 
Cpe (OGs Oar 05) e aituath bacea ie eek See ep 
SINE: 
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We may imagine that this curvilinear coordinate system moves con- 
tinuously + with the body as we pass from the original state By at time 
t—0 to the state B at time t. It will therefore form a curvilinear system 
in B so that, from (2.5) and (2.9), 

Y = (05, Dos Og, b)) . ae me ee eee 
Hence, the curvilinear system 6; may be identified with any convenient 


coordinate system in either By or B. 
We may define in B, contravariant components d6@' of a vector from the . 


relation (2.9). Thus 


i= SF dar dah Fas Se ae eee 
From (2.7) and (2.11) we have 
_ 06, OY esp: 
1oi—= — dy* Eee EL yay 
( ices dy oie . a ee eas 


so that d@’ remain the contravariant components of a vector under all 
transformations, including transformations from coordinate systems in 
B, to those in B. 

Using (2.9) and (2.10) the position vectors (2.1) and (2.3) of the points 
P, and P respectively take the forms 


r=r(6,, 05, 43), ] 
1» Pa» Os) (2.13) 
R=R(6,, 65, 63, t), 
and the displacement vector (2.4) becomes 
V=VW(07,05, Og.¢).02 >. Tae ee 
The line elements (2.2) and (2.8) which correspond to the vectors (2.13) are 
as =Gip dot dok, . . . . . . . (2.15) 
ds?=Gj, de? dd*, . . . . . . . (2.16) 


where g;;,, and Gj, are the covariant components of the metric tensors of 
B, and B respectively. The quantities G,,, are found for a given time ¢ 
and are functions of ¢ as well as 6,, 45, 43. 

The covariant components g,;, of the metric tensor of By, together 
with the corresponding contravariant components g‘*, are derived from 
the Euclidean space defined by the coordinates x;. These tensor 
components satisfy the appropriate laws of transformations for all changes 
of coordinates of the usual type, including changes to coordinate systems 
in B. Similarly, the covariant and contravariant components Gj, G‘* of 
the metric tensor of B are derived from the Euclidean space defined by the 
coordinates y,;, and they satisfy the appropriate laws of transformations 
for all changes of coordinates, including changes to coordinate systems 
in By. 


a a ae ge ee 
+ The values of 6; which define a point P, of By remain fixed with P, as it 
moves from its position in B, to its position P in B. : 
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Since both B, and B are situated in Euclidean spaces it follows that 
the components of the Riemann—Christoffel tensors for each body vanish, 
so that ' 


oR ijnr=9, Rip, Aa! a er Se ci) 


the prefix 9 referring to the unstrained body By. The tensor components 
oh j;,; are formed from the metric tensor components J,,,, 9” and their 
derivatives with respect to 6,. Similarly, the tensor components Riger 
are formed from G,,,,, and G”” and their derivatives with respect to 6,. 


§ 3. STRAIN TENSOR. 


The state of strain of the body B is characterized + by the difference 
of the squared line elements (2.15) and (2.16). We write 


ds*—ds?=2y,, a0! d6*, ee ee OSL) 
where 
VY wr=4(Gi,—9 sx)- 5 C 3 . . 5 (3.2) 


The quantities y;,, are the symmetrical covariant components of a tensor 
- which is called the strain tensor t. When y,,—0 for all points, ds =dsy, and 
the body is rigid. 

Since y,, are symmetrical it is always possible to find a coordinate 
system in B such that y;,—O(i=|=k). Since the coordinate system in B is 
then orthogonal, G,,=0 (¢=|=k), and hence, from (3.2), g;,=0 (t=|=4), so 
that the original directions of this coordinate system in By are also 
orthogonal. The directions of the coordinate curves in B which are 
determined by this coordinate system are called the principal directions 
of strain; and the strains in the principal directions, the covariant com- 
ponents of the canonical form of the strain tensor, are termed principal 
strains. 

With the help of (3.2) we can express the covariant components of the 
metric tensor of the body B in terms of the covariant components of the 
metric tensor of the original body B, and the covariant components 
of the strain tensor. Thus 


CoG ye te fe en ee a (3.3) 
and we obtain the determinant G of the elements Gj, in the form 


1__ 1 ,ijkplmn mx 
G= Fe” Se G1 GjmG en 


Pan ie ze wi Ee MG) UD imY ke n +GjmD enV it FI iG enVim 
+29 iV jmV en +I im iV em +9 enV Vim) FAY t1V3mV end bi 5 (3.4) 


+ The strain may be characterized in other ways as is done, for example, 
by Richter, H. (1949). 

+The strain in the neighbourhood of a point can be interpreted in a geo- 
metrical way by a strain “ quadric ’’, or more specially by a “ strain ellipsoid ” 
or “reciprocal strain ellipsoid’. See, for example, A. F. H. Love (1927), 
C. B. Biezeno and K. Grammel (1939). 
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where g is the determinant of the elements g , and where e* are completely 
skew-symmetric components of the e-system tensor, which take the 
values 0, +1, —1. 

In order to calculate the strain components in terms of the displacement 
vector we write (2.4) in the form 


R=r-+v—c > detiiadaal Raia 
and observe that 
E,=e;+V/;; oy AMIR Sts, eon a ae ET 
where 
E;=R ,, e;=r ;, Sidlaa Ge ay ee Ca 


e,, E; being the covariant base vectors in the original and strained bodies 
respectively. The comma denotes partial differentiation with respect to 
@;. Hence, since 
IR; - ey, Gy=E, . E,,; . . . . . (3.8) . 

we have, using (3.6), 

Gin=9in tes. Vinten - Vit Vii- Vix, 
and therefore 

Yn=s(e; e Vietex ° V cia . Vx) . ° ° . (3.9) 


The displacement vector may be expressed in terms of the contravariant 
base vectors e” of By. Thus 


v=»,,e" qt: topic, Bante Loe 
and 


V =| ie", SOUL ND ie OR Seek 


where the vertical line denotes covariant differentiation with respect to- 
By and means 


y 
Onl =m i pee VU». Sint SEE ee (3.12) 


The Christoffel symbols es) in the last expression are to be calculated 
q 


for By from the components g,,,,, g’””" of the metric tensor of B, and are 
given by 


ir ; 
— yrs y 
[s, Mi]yp=3HIsi,mo Ism, i Imi, s}- 


Introducing (3.10) into (3.9) the strain components (3.9) become 


(3.13) 


Vin=3( lp Hg] +e] M,ly)s Sots etre. es wits Mall te 
if we use the formule 


©7 CROP IN ae AO ae STS) 
where the Kronecker delta 6f has the usual meaning. 
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The strain components in (3.14) are expressed in terms of components 
of the displacement vector with respect to Bo and their covariant 
derivatives with respect to By). We may, of course, also express the 
‘strain-components in terms of components of the displacement vector 
with respect to B. For this purpose we write 

r=R—vte 
and we obtain the components of the metric tensor of By in the form 
ie Gie—E PV Ep. Vg tV it Ven 
and hence 
Vie t(E;. Vi, +E, . Vyj—ViG- Via + ws ++ (3.16) 
The displacement vector, in terms of the contravariant base vectors 
E’ of the body B, is 
Ve Ee eg ce nee, ah A cauk. ave Oc Lte) 
and its derivative is 
vo =VallE”. MER EEE Paes eae ee 
‘The double line denotes covariant differentiation with respect to B, that 
is with respect to.6; and G,,, so that 


T r r ian 
\ mY ee {nif ni, OF eS ete (3.1 9) 


where the Christoffel symbols (ee are calculated for B from G,,,,, G””. 
“Thus 


r 5 } 
{ }=Orts, mi], 
ae reed wee) 
[s, mi}=3{Gs;, piste kon ae Ae 


Since 
Eo Eh a6) oes boost 2s (3.21) 
-we find from (3.16) and (3.18) that the strain components now become 
Vin=H(Villet+Velli-V'U Velie) + 2 + + + (3-22) 
‘which gives the strain-components in terms of the components of the 
displacement with respect to B. 

The components of strain in terms of the displacement referred to any 
convenient set of coordinates in either B, or B may at once be written 
down from (3.14) and (3.22) respectively. For example, the strain- 
components referred to the cartesian system x; in By may be found from 
(3.14) by putting 0;=%;, and are 


ao, 14(B0)* 5 (@u)* 4 (22) 
ra Fe tsar) a ar) Vaz.) f? 


(52 +522) : (2 OV, | Wg AV, oS) 


¥23~ O\ 00, | Ot, T 3 0%, 0x, ' O%,0", | Ox, 0X5)" 


a 
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etc. Similarly, when at time ¢ the 6,-curves coincide with the fixed! 
cartesian y,-curves, so that 6,==y;, the strain components referred to ¥; 
may be found from (3.22). Thus 


EV (20\" (5) (5) | 
eee) an) en) J? 

(3 se) —3(Sc on ONION: a) 
D\ Gy, | -OYa/ -2\0Y, OY, Oa dYs O¥sOYs 


YON 
etc. 

The mixed components of the strain tensor can be formed in two ways. 
according to the choice of metric. Using the strained metric components. 
G'* we have 

eG STN SS Seer clea 
and using the components g‘* for the unstrained body, 


yeasty. o ae 2) SEE Oal eis 

Similarly, the contravariant components of the strain tensor may be 

formed in three ways. For our purpose the most convenient definition 
to take is 

thee * f= Gh — 9" Gt as. ee ee ee) 


Using (3.2) this becomes 


ek¥=1g"Gk\G,1—9 1) Be ety Bees l 
which reduces to 
ek—I(gik Git) ee AP en we 
or 
! Gthaa tha eh ee ee ee ee ee ae 
if we use the results 


9° j= 8, GHG. OF. UP gee eee 


§ 4. EQUATIONS OF COMPATIBILITY. 

The components y,, of the strain tensor, as found in the preceding 
section, are not independent but have to satisfy the condition that the: 
components of the Riemann—Christoffel tensors 9R,;;,; and R,;,, of the 
original and strained bodies must vanish. We may write the components. 
of the Riemann—Christoffel tensor of the strained body in the form 


i hee 
Rijn 2 (Gir gut Gyn, ir Gin, 37—Gy1, ix) 


+Gr{[m, li][n, kj]—[m, tk][n, lj]}. . (4.1) 
Using (3.3) the Cristoffel-symbols [¢, iy of the strained medium become 
[?, rs] =3(G Vir, ont Os: Arg, j=[, rS\ot+2y{;, rs) » + (4.2) 


where we have put 


Vi, rsj=2(Vir, or hes. r ea) 5 2 ‘ . (4.3) 
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If we now substitute (3.28), (3.3) and (4.2) into (4.1) we obtain 
Rjxi=oB geet it,gn +i, i Viel ata 

—2e"™([m, li]o[n, kj ]y—[m, tk]o[n, lj 1p) 

eg | [n, Kj lovim, wy +l, Woy tn, kj] 

—[%, Yor tm, sey —L ov tn, ws 
—4emnsj he : ) 
4e""tin, hj Loin, ay lm, liloytn, ag —l%, Y lo tm, ik] 
—[m, 2 Kovtn, ut 27 tm, KY [n, ki — 2m, ik 1Y En, ui 
+49" (4 KI (n, Bil — tm, sAY in, Bi) ee ER ay ML eS | 


Fig. 2. 


The conditions 
okt iki: . Rigi Q, - . : . . . (4.5)) 


used in (4.4) give relations between the strain-components which are 
called the equations of compatibility. As we know that there are only 
six independent components of the Riemann—Christoffel tensor, (4.4) and 
(4.5) reduce to six independent equations. 


§ 5. EXTERNAL Force. STRESS TENSOR. 


We assume that the body By is moved to its strained position B by 
the action of external forces which may vary with the time. At any 
point P of the strained body B we construct a parallelepiped which is 
bounded by the faces 0;=constant, 0;+d0'=constant, and a tetrahedron 
whose edges are formed by the coordinate curves PP of lengths ds;, and 
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the curves P,P,, P,P3, P3P, (see fig. 2), which in the limit are defined by 
the three vectors 
ds;—ds,=E, d6/—E, d0! (i, jnot summed). . . . (5.1) 
The surfaces 0;—constant of the tetrahedron have areas 
1d8,=4 V/(GG")de) de® (i= J =|= &, ¢ not summed), . (5.2) 
and the area of P, P,P, is denoted by 3d8. 
If n is a unit vector normal to P,; P,P; then 
n dS=(ds,—ds,) x (ds,—ds,), 
and, using (5.1), this becomes 
ndS=(E, x E,,)d0/ d6* 
") it (j, b=; 252,83: 3, 10 eee 
=e ;;,E!' db! do* 
where 
€;j,E'=E, XE, . ry . . e ; . (5.4) 
€jp—= VG if i, 7, k is an even permutation of 
the numbers 1, 2, 3, 
=—V/G if 7, 9, k is an odd permutation of (5.5) 
the numbers 1, 2, 3, x ae % 
==() if any two of the indices 7, 7, k are 
equal and not summed. 


Hence from (5.3) and (5.2) 


Ei 
Yay Ge 

If n; are the covariant components of n, with respect to the base 
vectors E’, it follows that + 


Ri/(Gi)\dS=dS &.-<ae. Se 


A stress vector t is associated with the surface P,P,P, of the tetra- 
hedron PP,P,P,. The stress associated with the surfaces 6,—constant of 
the tetrahedron can be represented in magnitude and direction by 
quantities —t;. The quantities —t; are not invariant under general 
transformations of coordinates. This fact is apparent from the notation 
‘since vector invariants do not contain suffixes. The suffixes, however, 
do not indicate that the quantity transforms according to the covariant 
type of transformation. 


By considering the equation of motion of the infinitesimal tetrahedron 
PP,P,P, we have, in the limit, 


tdS=t, ds ;. 5 7 4 4 e ° ° . (5.8) 

a ee a a et ea 

+ We remember that the summation convention is suspended when an index 
appears more than twice. 


S;. 


t 


(5.6) 
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Volume forces and mass-accelerations acting on the tetrahedron do not 
appear in this equation since they are of a higher order of smallness than 
the surface forces. 


From (5.7) and (5.8) we find that 
t= Ln t/G". ; . 5 . 5 ° 5 . (5.9) 
i 


Since, under general transformations of coordinates, t is an invariant and 
n, are the covariant components of a vector, it follows from (5.9) that 
t,/G" transforms according to a contravariant type of transformation. 
We may therefore write 

ty Glee "Er Ens Oe eg 6.10) 
where 7’ and 7, are the contravariant and mixed components of a tensor 
of the second order which is called the stress tensor. Later we will prove 
that the contravariant components are symmetric so that we are allowed 
to write the mixed components as 7}, and not 7*,, or 7j!. 

The stress quantities t; exerted by the body in the region 0;>0 on the 
body in the region 0;<0, across the 6;-surface, has positive values for its 
components when 7’* are positive. 

The covariant components of the stress tensor may be defined by the 
equations 

Ty ean int Se een ee ee ete ee, (Oe 11) 


From (5.10) we see that the contravariant components 7‘* of the 
stress tensor and the mixed components 7’, of the stress tensor are related 
to the quantities t,,/G". The covariant components 7,, of the stress 
tensor cannot be related to the quantities t;,/G" in a simple way and are 
therefore, in general, of less importance. 

From (5.9) and (5.10) we obtain 


teary Betis ee Bea te. (5212) 
so that if 
6S at ee eee cad tgs (0213) 
then 
ph ila) po ipa ee ht a) ney 27. (0.14) 


and ¢*, ¢, are the contravariant and covariant components respectively 


of the stress vector t. 
The three stress quantities t; in (5.10) may be put in the form 


= Give) ixn_¥s 6.15 
az (ZB): Ree . . . . . ( . ) 


where E,/\/G,, are unit vectors along the coordinate curves. The 
physical stress quantities are therefore 


i (Gi) = ee ee nee teen (D6) 


These quantities are not, of course, the components of a tensor. 
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§ 6. Equations oF Morrion AND EQUILIBRIUM. 


We consider the infinitesimal parallelepiped in fig. 3 whose faces have 
areas dS, and we consider the equations of motion of all the forces acting 
on this elementary volume. The forces on the faces @;=constant through 
P are, to the first order, —t,; dS; (¢ not summed), and since 


dS = /(GG")dei dé” (¢ not summed, i=[=j=|-k), . (6.1) 
we may write this as 
—T,d6id® (i= jE), 
where 
T.=/(GG")t,—+/(G)r#E,=1/(G)AE*.. . tes 
The forces on the faces 6;=constant and on the faces 6;+d@'=constant 


are shown in fig. 3. 


Fig. 3. 


(Ts +Ts,3 20°) 20267 


(h+T,.de)a9'de* 


1,106 aed? 


. The body forces acting on the elementary parallelepiped whose volume- 
is dr are pF dr or pQ dé dé? dé3, where 


Q=FV/G, +, sdye wks saa en 


and where p is the density of the strained body B. 

We now take moments about the centre of gravity of the elementary 
volume, and, since in the limit we may retain only first order terms, this 
gives 
=) 


E,xT,=0. . . . ° . e . . (6.4) 
From (6.2), (5.4) and (6.4) we obtain 

€,,17 *E'=0. 
It follows immediately from this equation that , 


ik pki diy et oul ieee ae 
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so that the contravariant components (and therefore also the covariant 
components) of the stress tensor are symmetric. 

Since the components 7‘* are symmetric it is always possible to find a 
coordinate system in B in which 7‘*=0 (i=|=k). The directions of the 
coordinate curves in such a system are called the principal axes of stress 
and the stresses which correspond to these axes are called the principal 
-StLeSSES. 

The remaining equation of motion states that the (vector) sum of all the 
forces acting on the elementary volume equals the mass-acceleration of 
the volume, where the acceleration may be taken to be that of the centre 
of gravity. Hence, in the limit, (see fig. 3) 

T;, -+pQ=ph, et dees 3) i ae OO) 
where 
h=f,/G, ee Ue eee ce oe (On) 
and f is the acceleration vector. 
When the body is in equilibrium f=0 and 


T;, + pQ=0. 5 . . . . . . (6.8) 
We may express the equation of motion in alternative forms. If 
Fab =F ,E : 
| (6.9) 
f=fE=fE", J 
then 
= iE ae Ee 
ap he (6.10) 
h=/E=1,E% 
where 
i i G, =, G, 4) 
. He oh vy Aut E5 (G20) 
ia fta/ Ge. W=f A/G. 
Also, if hints 
otk—7tk,/G, o=7,7 G, Tore a ae Onde) 
then, from (6.2), 
To*E =oEF 9. . .i. 9°. (6.13) 


Tt should be noticed that the quantities defined in (6.11) and (6.12) are 
not components of tensorst. Using (6.9) to (6.13) the equation of motion 
(6.6) becomes 


a i {i fo*-beQ*= ph 6 wey ee A (6.14) 
‘or 


ai. {sport ei eh Sere oe et (0515) 


Asm bach Ace 
\a}=Ve Co 


+ They are sometimes called relative or weighted tensors. 


‘Since 
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these equations may be reduced respectively to 


rik - rit 4 Lat ttt eRe OE AGRE TOY 
hin dapat{it Tit pFp=pfy + + + + (617) 

These may be written more compactly in the tensor forms 
rthed | Coot S Vert.) LA ee meee i = ee TREES ES 


7t||,+pF,.= phy Ss inate th sy fe 8 


where the double line again denotes covariant differentiation with respect 
to the strained body B. 

With the help of (5.11) it is possible, of course, to write the equations 
of motion in terms of the covariant components of the stress tensor but 
this form appears to be of little interest and will not be given here. 

The equations (6.6) and (6.14) to (6.19) are alternative forms of the 
equations of motion. The most convenient form to use depends on the 
type of problem which is being examined but the equations (6.6), (6.18) 
and (6.19) are particularly simple in form. 

Finally, the conditions at the boundary surface of the body at which 
the surface forces are prescribed require that, 


t=?” 5). ee 


where P is the surface force vector. If we express the surface forces P in 
terms of their components, 


Ba PkE = PE®. 4 pee 
and use (5.12), the boundary conditions have the alternative forms 
| chp Pk -9 ois. Se at ee ae ana e 
or 
Ane=P ip wl Jan) Sawicn epee 


§ 7. Virtual WorRK. 


When the body is in the strained state B at time ¢ the displacement v 
of each point P which was originally at P, is 


v=Vv(6;, 09, Og, ¢). 


We now consider any other geometrically possible motion which is such. 
that, during the same time ¢, the points P, move to P’ and form a strained 
body B’, such that the displacement of P, is of the form 


v’=v(6,, 05, O5, t)+-ew(01, 95, 95, t), 
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where w is an arbitrary vector and ¢ isa real number which we shall take 
to be small, so that its squares and higher powers can be neglected compared. 
with «. We call ew a variation or virtual displacement of v and write 


ov=v' —v=ew. 


This corresponds to the usual, definition of variation in the calculus of’ 
variations}. | 
The variation of v,; is defined to be 


d(V.j)=V" 


iV, =ew, ;=(ew) ;=(6v) ;, Ao ee 


so that 6 and partial derivatives with respect to 6; are interchangeable. 
The variations of the covariant base vectors E; are given by 


dE ;=6(e;+-v_;)=(8v) |; 


if we remember that e; does not depend onv. The variations of the scalar- 
products of the base vectors are found to be 


Hence, using (7.1) and (3.2), the variation of the components of the strain 
tensor are 


By ,,=—6G,,=E; . (dv) ,+-E, . (dv) |; . . . . (7.2) 


since 6g ;,=0. | 

We now consider the whole strained body and calculate the work 
done, during a virtual displacement, by the body forces and the surface 
forces acting on its boundary. This is called the virtual work and is. 


| 


defined fT to be es 
seA=| P.dvdS+ | BEL Sve ee het) (E3) 
SS) v 


where dS is an element of the bounding surface and dz a volume element.. 
In (7.3) the surface integral is to be taken over the whole of the bounding 
surface S. The volume integral is throughout the volume of B. In order 
to be able to evaluate the surface integral in (7.3) we assume that the 
surface forces are prescribed over parts of the surface and that the surface 
displacements are prescribed over the remainder of the surface. If the 
variation of v is then chosen to be zero where the surface displacements. 
are prescribed, the surface integral is zero over these parts and may be 
evaluated over the remainder of the surface over which the surface forces. 
are known. 

Using the boundary conditions (6.20) the expression for 6*A may be 
written 


stA— | t.dvdS+ | PERO VT ee! (7.45 
s : 


+ See, for example, Courant and Hilbert (1931). , 
+ The star on 5* indicates that 5*A is not the variation of a function A. 
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From (5.12) and (6.2) 
n 


ty Gan . . . . . . . . (7. 


~l 
Or 
— 


so that (7.4) becomes 


T 


Bea =| Gon bv dS + | oF . 8v dr. 


‘The surface integral may be transformed by using Green’s theorem in the 
form 


PABVG ie 
| Brn, dS = | AB VG Siege (7.6) 
8 oar V G 
- where B’ are the contravariant components of a vector. Thus 
= ' T. 7 5vV 27 r fa dr 
S*A = ik eae dr + " pF. bv dr=| {(T,+eQ)v+T;. (8¥),3} 7G: 


and using (6.2) and (6.6) this reduces to 


S*A = | {pf. v+7%E, . (Sv), fdr. DA 


Lf if 


Since 7’* are symmetric components of a tensor 


S*A = | {pf. 5v+47"[E, . (8v),+E;. (Sv), .]} dr 


= | | of. dv+7ydvG do ded®, . . . (7.8) 


if we use (7.2). 

Subject to boundary conditions at the surface of the body, the virtual 
displacement is an arbitrary geometrically possible displacement, so that 
equation (7.8), as is seen from its derivation, can replace the equations of 
motion of the body. 

When the body is in equilibrium, f=0, and (7.8) reduces to T 


ee J { rik(Sy.)4/G do! de2d63. =. . . . (7.9) 


‘This represents the necessary (and sufficient) conditions that the body 
is in equilibrium and may replace the equations of equilibrium. 
For a rigid virtual displacement dy ,;,,=0 and therefore 


S*AS0.0, 2 soa nk se ee DE 


This states that for a body in equilibrium the virtual work of the external 
forces is zero for any rigid virtual displacements. 


+ This result is a general form of the well-known principle of virtual work of 
elastic bodies. 
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§ 8. SrrAIn-ENERGY FUNCTION. 


We now confine our attention to a body By which at time t=0 is 
unstrained and unstressed +, and we assume that B is obtained from B, 
entirely by the action of external forces t. If B is able to recover its 
original size and shape completely when the external forces cease to act 
the body is said to be elastic. This implies that when all the components 
of the strain tensor vanish then all the components of the stress tensor 
vanish and vice versa. 


An element of mass of the strained body B at time ¢ is given by 
dm=p»/G de dé dé, Segtinyit baad to WRC SEL) 


and the principle of conservation of mass in a virtual displacement is 
expressed by 


d(dm)=0. {Ray eee ee ee Coe} 

We now consider the actual motion of the body B which takes place 
“Lae : ¢ 0 : 

in time dt so that the displacement vector increases by dv= s , 9; being 


fixed. Let dK denote the increment of kinetic energy of B during the 
time dt. The kinetic energy K is given by 


dv\2 
c= 1 = ae ¢ 
K=1{ (3) Cie terre ahi d  Atated hin id90( 8.9) 


Ov Ov 
Te | 5d (=) ay 


’ Ov\ . ‘ 
since the mass element dm is unaltered. But d (3) is the increase of 


so that 


velocity in time dt and therefore 


Ov 0?v 


Hence 


aK = | £.dv dn= | || (f. dv)p/G db de? aes. BW (4) 
B 


The actual increase in the work done d*A by the external forces acting 
on B, during time dt, may be obtained by putting v=dv and 


ay; 
by i, = dy ik— “i dt, 


BS pO EE IE met PB Fae Ea 2 ea lo 
+ This means, for example, that no stresses are induced in By by its 
manufacture or by the action of body forces. 
+ Strains produced without external forces but by changes of temperature 
are not considered. 
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in (7.8), where dy,,, are the increases in the components of strain in time dt. 
Thus 


d*A= { i | (pf. dvr dy,),/G dO! de2d63, . . . (8.5) 


The fundamental energy-conservation laws of thermodynamics for 
reversible processes may be written 


{ T do dm-+d*A= | dU dm-+-dK, =. i. (8.6) 
B B 


where T is the temperature of the mass element dm, do the increment of 
entropy per unit mass and dU the increment of internal energy per unit 
mass. By using (8.4) and (8.5) this equation reduces to 


[ff ao doyp ye aon doe a= fff =a) VG dO% a62 abs. (8.7) 


Since we may consider any arbitrary volume, it follows that 
(@U—T do)p=—7* dy} 5 uae 2b eee ad 


if we make suitable assumptions about the continuity of the integrands. 
If the changes which have taken place in time dt are reversible and 
isothermal, then T is constant, and 


pile 
dW==d(U—To)=— dy - - - . «= (8.9) 


where dW is a complete differential. The quantities dy, are also, for 
fixed values of 6;, complete differentials. Since we are assuming that the 
body is elastic, the changes which take place are reversible, and the 
condition of the body B is independent of the actual path taken to reach 
the state B from By. The differentials dy;, in (8.9) may therefore be 
considered as arbitrary so that, since dW is a complete differential, 


min oe re, ere ee FN 
Ovin: 


where W is a function of y;, and all changes of state are isothermal. 
Similarly, if the changes in an elastic body are adiabatic, so that 
do=0, we find from (8.8) 


he ap ek Sy teres Ye 


where U is a function of y;,. 

The function W and U are called strain-energy functions and play 
the réle of potential energy per unit mass for the isothermal and adiabatic 
cases respectively. 


+ Under these conditions dU and de are complete differentials. 
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The relations (8.10) and (8.11) between the contravariant components 
of the stress tensor and the strain-energy function is the most convenient 
for our purpose, but it is of interest to derive an alternative form t. 
The mixed components ¢, of the strain tensor are defined in (3.23) and. 
from this definition 

Vat 1 R= 1G w+ 2Y i) 
Hence 
dy 7,5=Gyj dey + 2¢ dy; 
and therefore 
GMO Fe ey eG GO, dade... . (8.12) 
From (8.10) 
paW=r"* dy ,—Gmrh, dy pi, 
and with (8.12) this becomes 
p(di—2e!) dW=7*, de”, 


so that [, for the reversible isothermal case, 
: oW ;oW see 
=e (Sm 245), . . . . ° . (8.13) 
with a similar result with U written for W when the changes are reversible 
and adiabatic. 


§ 9. VARIATIONAL, WORK AND ENERGY PRINCIPLES. 


In this, and in subsequent sections, we shall consider an elastic body 
which undergoes isothermal changes of state and use the strain-energy 
function W. The corresponding results for adiabatic changes may be 
obtained in a similar manner. 

If W’ is the total strain-energy of the whole body, then 


w'= | [[pWwyG ae ae aoe, eet 25 (9:1) 
and, remembering § (8.2), 


sW’= { i { (SW)p/G dot dé? a9? 


es { { i 5 Oyule VG do! di de, hoe (9.2) 


+ Other relations between the components of the stress tensor and the 
strain energy may be found by using the expression (7.7) for 6*A instead. 
of (7.8). 

+ If our coordinates 6; are identified with a general coordinate system in the 
strained body B at a given time ¢, then (8.13) agrees with a similar result 
obtained by a somewhat different argument by F. D. Murnaghan (1937). On the 
other hand, if we take 0; to be a general coordinate system in the original body Bo 
the result (8.10) is the same as that given by Murnaghan when his two coordinate 
systems coincide. : oe 

§ The 8 used here is the usual variation symbol of the Calculus of Variations. 


2B2 
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Hence, from (7.8), (8.10) and (9.2), we have 
S*A—SW' {lf (f. Sv)p\/G dei de? de®. . . . (9.3) 
We integrate this last result between two arbitary times ¢) and ¢,. Thus 
[ora di— ["Sw’ at+ fi {ff (f. 5v)p\/G dol de? de dt. . (9.4) 
. to to 


The last integral may be transformed by partial integration. Using 
2: f : 
f= a and remembering that the mass element p./G dé" dé? dé? is 


constant, we obtain 
hy 


j J J J (f-8v)pyGae dt de di— | { { | (3. iv) p\/G do! dé? aoe | 
to 
“pS Movecananaea 


We now suppose that the variation 6v is such that it vanishes at time f, 
and t, so that 


ty ty 
| | I i (£. 8v)p+/G do! de? dé? dt=—[(8Kdt, . . (9.5) 
to 


to 
where K, the kinetic energy of the strained body B, is given by (8.3), and 


= [5 -3(Gp) eve ae aoe ae ig a Le 


Equation (9.4) now becomes 


ty ty ty 
S*A ‘dias (@8W’ di—2.{iodKidis ply aan 


to to to 
If the external forces can be derived from a potential energy function 
Q then 
d*A=—S8Q, a eee ek re eS 


and (9.7) reduces to the simple form 


th ty 
8(K—W’— 2)dt=8 { (Ki— Wi=,Q)di=0, . eee) 
to to 

where K—W’— 2 is the Lagrangian function. 

Equation (9.9) is the analytical expression of Hamilton’s Principle, 
which says that the variation of the integral of the Lagrangian function, 
between times f) and ¢,, takes a stationary value, provided that the 
variation of the displacement vector is taken in such a way that it vanishes 
at the times f, ¢,. 

If the strained body B is in equilibrium then f=0 and the virtual work 
(9.3) of the external forces acting on the body becomes 


d*A=dW’, 
and, if (9.8) still holds, 
SOW! iQ) 0. njceeeahney copie daar k Gs10) 
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This states that the variation of the total potential energy W’+ 2 has a 
stationary value. . 

The principles (9.9) and (9.10) were found from the expression (9.3) 
for the virtual work of the external forces. We may also derive from 
(9.3) expressions for the actual work performed in time dt by replacing 
dv by the actual displacement dv. The virtual work of the external 
forces 6*A now becomes the actual increase d*A of work performed in 
time dt, and 6W’ becomes dW’, and therefore (9.3) is now 


d*A—dW' + { | i (f. dv)p\/G 01 ded. . . . (9.11) 


If we use (8.4) this reduces to 
9 


Read Wied Keene hs ee me nO DY 


This states that the work done by the external forces in time dt is equal 
to the sum of the increase of the kinetic energy and the increase of the 
strain-energy. 


§10. SrrEess-STRALIN RELATIONS. 


We have seen in §8 that, under certain conditions, the stresses which 
occur in an elastic body are related to the function W by the equations 


ano W) 
WE a ay 10.1) 
; oye: 
_ and 
Wie | | | BV OdOLdhsdge me oe 4. (10.2) 


represents the total strain-energy of the strained body. Since mass 
elements are conserved under straining, we have 


PoVI=pvG;, eae eis A) eae e's (10.3) 


where py is the density of the unstrained body. Hence equations (10.1) 
may be written 
a/(G)7**=pyv/(9) oe ee eos) (1054) 
Vik 
At this point it is necessary to make further assumptions if we wish to 
obtain an explicit relation between stress and strain. The validity of the 
assumptions when applied to physical bodies must be tested by referring 
toexperiment. The relation between stress and strain, in the theory which 
has been developed so far, is of such a form that a strain-energy function. 
W exists with the property (10.4). We therefore assume a suitable form 
for W and derive the stress-strain relations from (10.4). 
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Since the density p) of the unstrained body, and W, are scalars, the 
quantity pyW is invariant under all transformations of coordinates and 
we assume that W can be expanded as a power series ¢ in the strain 
components y;,. Thus . 

poW =E+ Bihy + Bry yp tReet ymnt---» (10.5) 
The functions E**, F’*”, ... do not depend on y,;,and since W is invariant 
they must be components of tensors. If W is taken to be zero when the 
body is in the initial state B in which y;,—0 then E=0. 

From (10.4) and (10.5) we have { 


J (Gig)tt= Ek Bity Beem y yyy. + + . (10:6) 
where we have put 
Kikr— Pike Prik, 
(1087) 
K tkamn — iklrmn | Blrikmn | Fmnlrik | 


Since we have assumed that our material is elastic there is no stress 
when the strains are zero provided there is no initial strain. Hence 


H*k—0. 
The expression (10.7) gives us a general relation between stress and 
strain, and we call the coefficients E’*", Kit", . . . the elastic coefficients 


of the body. They depend on the properties of the body in the unstrained 
state. 

The quantities EK’ are coefficients in an expansion of the strain- 
energy function W in terms of the components of strain and therefore 
refer to an elastic body which undergoes isothermal changes of state. 
When the changes are adiabatic the strain-energy function is U and the _ 
corresponding elastic coefficients are, in general, different, but for many 
purposes these differences are not great and can be ignored. 

From the symmetric properties of the components of the stress and 
strain tensors, and the relations (10.7), we deduce the following symmetric 
properties : 


Bik Fo kilr — Byikri Erik ee os. MO (10.8) 
with similar properties for Kim, , 
Because of the symmetric properties of the elastic coefficients, we see, 


from (10.4) and (10.6), that the strain-energy (10.5) may be expressed in 
the more convenient form 


PoW = SE ary yt PE ey LY mn * a (10.9) 

In order to examine the properties of the elastic coefficients in more 
detail we refer all quantities to the rectangular coordinates x, in the 
unstrained body By so that 6;=2,; and g=1. In these coordinates the 
tensor components E‘*ir, Hikirmn |__| are denoted by ci*ir, cikirmn, 
and the covariant components y,,,, of the strain tensor by e The 


mn 


+ This, of course, is not the only possible assumption that can be made. 
+ In using (10.4) we suppose that y;, is distinct from yp. 
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contravariant components of the stress tensor are denoted by t’* and these 
are referred to coordinates in the strained body B which, in the unstrained 
body Bo, were the coordinates x;. The strain-energy function then becomes 


PoW =—2c'# "ee, 4c rnc elma ts +s ..-» ~ (10.10) 

and 
| EN (G)=ac'"te, citing et et pet ee (10.11) 
If the elastic coefficients c’*’", c’krm,... are independent of the 


orientation of the rectangular axes x; at each point we say that the body 
is elastically isotropic, or more simply, isotropic. Otherwise the body is 
said to be zolotropic or anisotropic. 

If the density of the original body By is constant and if the elastic 
coefficients ci*", c¢ikirmn, |. are constants, we say that the body is 
élastically homogeneous. In this case the covariant derivatives of the 
general elastic coefficients E'**, Hikrmn, ||. are zero, so that 


ee El Vee rey io vaca, (10212) 


In the remainder of this paper we shall restrict our attention to bodies 
which are homogeneous i in the above sense so that equations (10.12) hold. 
Similar results, however, hold for non-homogeneous bodies. Because of 
theoretical and experimental difficulties we restrict our attention to those 


materials for which the coefficients E‘*”"”, ... are sufficiently small, so 
that we can express the strain-energy W approximately in the form 
| PoW =3 Ey ay p= Ferre Cr,» tt ying (el OS) 
‘with the corresponding stress-strain relations 
eG ety), 2, fuvide  ane (10514) 
Ay (Giese te. Rat a oe leet cme! Meh) 
The strain energy in this case can be also put in the alternative form 
Bees eel. hn) cages (L0v1G) 


§11. Exastic COEFFICIENTS. 
In general, remembering the symmetric properties 


ctklr— ckilr — cikrl — clrik 2 ky se" CE ee ace erie) 


we have 21 elastic coefficients c’" for an sxolotropic body. If the 
coordinate system a; is changed to another rectangular cartesian system 
z; then the coefficients c'”” become d‘ and, in view of their tensor 


character, 
Oz; 0%, O28, Oz 
qitie — ee Cc cet nein ci LL 
Fagen Guay iit?) 
Symmetry of structure of the general eolotropic body will introduce 
relations among the elastic coefficients c‘”" which may be found directly 


by using (11.2) and suitable transformations of Soominrrrees) but we will 
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not examine these in detail here. When the body is homogeneous and. 
isotropic the coefficients c“” are invariant under all transformations to- 
other rectangular cartesian coordinates and it can be shown that the 
coefficients reduce to 


cmnst — ASMMS st + SMH ME SHMH"S, “4 ee (1 l .3) 
where A, » are constants. Hence + 
Hiklr— \gikg@ u(gilgktgitgk!), i Mil 8% ele ag (11.4) 
or ’ 
Bip glgh Lott fe et gig. ia eh oe 
where 
relat 
1-27) - 


From (10.14) and (11.5) we find that the stress-strain relation for an 
isotropic body takes the form 


ViG)r*=nvo){a"a 0% + a0" bY A ee Cae 

or 
Er se 2 AS < 
VG) nv goteaeyt + ah. 2, weabhbaa 


If the coordinate curves 0; are along the principal directions of strain 
then y jx, Fix, Ging, G'* are all zero when 7=|=k. Hence, when the stress— 
strain relations are of the form (11.6), 7*=0 (i=|=k), so that the principal 
directions of stress in the strained body B coincide with the principal 
directions of strain. 


REFERENCES. 


Brezeno, C. B., and GRaMMEL, R., 1939, Technische Dynamik. Springer, Berlin. 

Britiouny, L., 1938, Les Tenseurs en Mécanique et en Elasticité, Masson, Paris, 
p. 226. 

Courant, R., and Hitpert, D., 1931, Methoden der Mathematischen Physik. 
Vol. 1. Springer, Berlin. 

Devuxer, E. A., 1941, Deutsche Mathematik, 5, 546. 

GLEYZAL, A., 1949, Quart. App. Math., 6, 429. 

Lovn, A. E. H., 1927, A Treatise on the Mathematical Theory of Elasticity. 
Cambridge University Press. 

Mitne-THomson, L. M., 1948, Proc. 7th International Congress for Applied 
Mechanics, 38, 

Murnaauan, F. D., 1937, Amer. Jour. Math., 59, 233. 

Nevuser, H., 1943, Z. angew. Math. Mech., 23, 321. 

Ricuter, H., 1948, Z. angew. Math. Mech., 28, 205. 

Ricuter, H., 1949, Z. angew. Math. Mech., 29, 65. ; 

Trerrtz, E., 1928, Handbuch der Physik, Vol. 6. Springer, Berlin 1939. 


} This result may also be obtained by assuming that the strain-energy 
function W is a function of the strain invariants. See, for example, F. D. 
Murnaghan (1937). 


[33744 


XXVIII. The Emission of Short- range Alpha Particles from Light 
Elements under Proton Bombardment.—Il. Further Observations on 
the Reaction B(px)'Be. 


By W. E. Burcuam* and Joan M. Freeman f. 


[Received December 23, 1949.] 


SUMMARY. 


The energy release in the reaction B(pa)’Be has been measured, using 
targets of known thickness, and found to be 1-147-+0-010 keV. Evidence 
for the existence of an excited state of energy about 400 keV. in the 7Be 
nucleus was obtained. The cross-section for the reaction was found by 
comparing the yield of alpha particles with the yield from the B(na)*Li 
reaction observed with the same target. A value of 4:5 10-27 cm.2 at 
330 keV. proton energy was found. 


§ 1. INTRODUCTION. 


IN a previous paper (Burcham and Freeman 1949, subsequently wierd 
to as I.) measurements of the energy release and cross-section for the 
reaction 

WB te Be tHe. Nes ci ose es (1) 


_ were reported. An improved cross-section determination has now been 
made by comparing the yield of reaction (1) with that of the slow neutron 
reaction 

WB+ In>"*Li-+*He. imam oor rec (2) 


This work, together with further observations on the energy release 
and excitation function for reaction (1), is described in the present paper. 


§ 2. ENERGY RELEASE IN THE REACTION !B(px)‘Be. 
(a) Q-value for the ground-state transition. 

The energy release reported in 1. was 1-11+0-06 MeV., while in a 
similar investigation Chao, Lauritsen and Tollestrup (1949) found a value 
of 1:148-+-0-005 MeV. The large error in the former determination arose 
partly from uncertainty in the thickness of the isotopic targets used, and 
in order to eliminate this error a special study of the alpha-particle group 
emitted from’ boron targets of different thickness has now been made. 
A set of targets each of area 0-20 cm.?, and of thickness known from the 
ion current in the electromagnetic separator, was prepared ; the density 
of the deposits ranged from 30 to 300 micrograms per cm.?. These 
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targets were bombarded with 432 keV. protons, and the alpha particles 
emitted at 123° to the proton beam were observed in the magnetic 
analyser as described in I. Some of the curves thus obtained are shown 
in fig. 1 ; the energy scale for these curves was obtained from a calibration 
of the magnetic analyser carried out by means of protons of energy 
known to +5 keV., scattered from a thin gold film as described in. I. 
In the experiments with boron a fluorescent screen of circular aperture 


Fig. 1. 
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Alpha-particle energy (keV.). 


Alpha-particles from the reaction 1B(pa)’Be. Proton energy 432 keV. ; 
angle of observation 123° ; target thicknesses 43, 63, 170, 310 and 500 keV. 
for alpha particles. 


1-0 cm. in diameter was used, and the resolving power of the analyser, 
as defined in I., was 25 (width 66 keV. for alpha particles of energy 
830 keV.). From the known analyser resolution and the excitation 
function for the reaction a curve representing the energy spectrum to 
be expected from a target of given thickness was constructed. For this 
purpose it was assumed that the alpha particles emerging from a particular 
depth in the target gave a spectrum which could be approximated by a 


Emission of Short-range Alpha Particles under Proton Bombardment 339 


triangle of height proportional to the cross-section for the reaction at the 
_ average proton energy for this depth in the target, and of width (at half- 
maximum) determined by the resolution of the analyser. A series of 
such triangles was drawn with 20 keV. intervals, the first representing 
alpha particles coming from an average depth of 10 keV., the second 
representing alpha particles from a depth of 30 keV., and so on, for a 
target of total thickness equivalent to 300 keV. (0-2 mgm./cm.2), as shown 
in fig. 2. An integrated curve which would be expected to represent the 
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900 


observed spectrum was thus obtained (full curve). Similar curves, shown 
dotted in fig. 2, were drawn for targets of equivalent thickness 40, 80 and 
' 120 keV. The intensity falls off on the left-hand side of the peak due to 
the effects of the exponential form of the excitation function and of the 
decreasing resolution width. It can be seen that for targets of equivalent 
thickness 120 keV. or more the peak of the spectrum occurs at an energy 
60 keV. less than that for an infinitely thin target, and that only the first 
120 keV. equivalent thickness of the target contributes to the peak. 
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The position and intensity of the latter should thus be independent of 
the thickness for such targets. For thinner targets the intensity would 
be expected to decrease and the peak to move to the right. A comparison 
of the constructed curves and the experimental curves of fig. 1 shows that 
the latter set, although somewhat narrower, exhibit the general features. 
expected, and justifies the procedure of deducing from the constructed 
curves the depth in the target in keV. equivalent for which the peak of 
the alpha-particle group from each of the experimental targets would 
be expected to occur. The stopping power of the boron targets for 
alpha particles of about 800 keV. energy was assumed to be 1450 keV. 
per mgm./cm.2, which is consistent with the data of Madsen and 
Venkateswarlu (1948) on the loss of energy of protons in beryllium*. 
The results are set out in Table I., which gives the actual target thickness. 
in mgm./cm.?, the effective thickness allowing for inclination @ of the 
target to the proton beam, the corresponding thickness in keV. equivalent 


TABLE I. 
Target thickness Depth of peak Intensity at | Energy at 
from surface | peak (observed) peak 
mgm. /om.? mgm./cm.” keV. (keV.) alpha-particles | (observed) 
sia! cos 6 equiv.| (theoretical) |perl0!protons} (keV.) 
0-300 0-345 500 60 5:7 805 
0-185 0-215 310 60 5-6 804 
0-165 0-190 270 60 4-9 814 
0-105 0-120 170 60 5:1 815 
0-038 0-044 63 32 5-0 819 
0-026 0-030 43 22 4-2 825 


for alpha particles, the depth in the target for which the peak would be 
expected to occur, and the alpha-particle intensity and energy observed 
at the peak of each group. 

When the alpha-particle energy is plotted against the theoretical depth 
of the peak from the surface (column 4) the result shown in fig. 3 is 
obtained. The extrapolation to zero target thickness gives an alpha- 
particle energy of 833 keV., with an uncertainty in the extrapolation 
of 3 keV. This figure must be further corrected to allow for the passage 
of the alpha particles through the thin layer of carbon deposited on the 
surface of the targets from pump-oil vapours. The thickness of this 
layer was found by determining the yield of radioactive N1°, formed in 
the C(py) reaction, from the boron targets in comparison with that 
obtained from a thick carbon target. The average value of the thickness 


* Note added in proof.—The stopping power of beryllium for protons of 
velocity equal to that of 800 keV. alpha particles can be obtained directly 
from the recent work of Warsaw (1949, Phys. Rev., 76, 1759) 
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found for the boron targets was equivalent to 3-L1 keV. for protons ; this 
correction has been made to the incident proton energy and a correction 
of 12+4 keV. to the emergent alpha-particle energy. With the 
resulting proton energy of 429+7 keV. and alpha-particle energy of 
84547 keV. at 123° the Q-value for the ground state transition in the 
B(pax) reaction is found to be 


Q=1-147-+0-010 keV., 


in good agreement with the work of Chao, Lauritsen and Tollestrup (1949). 
The error in this Q-value is less than that calculated directly from the 
stated errors in the alpha-particle and proton energies, since calibration 
errors affect these two energies in the same sense. 
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Reaction 1°B(pa)’Be. Extrapolation of observed alpha particle energy to zero 
target thickness. 


The curves of fig. 1 show that the shape of the high energy side of the 
alpha-particle distribution is independent of target thickness. An 
estimate of the Q-value from an “extrapolated ”’ energy derived by 
drawing a tangent to the observed spectrum gave 1-:157-+-0-020 MeV. ; 
this is less accurate than the preceding value owing to greater uncertainty 
in the analyser calibration for “* extrapolated ” ranges. 

(b) The first excited state of “Be. 

It was pointed out in I. that the similarity between the nuclei “Li and 
7Be suggested that there might be a low-lying state of “Be analogous to 
the 480 keV. state of 7Li. A search (I. fig. 10) for alpha particles emitted 
in the 1B(px)?Be reaction which would correspond to transitions to 
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such an excited state showed that the intensity of such a group was less. 
than 10 per cent of that of the main group for a proton energy of 230 keV. 
If, however, the a priori probability of emission of the two groups is the 
same, barrier penetrability would reduce the intensity of an alpha-particle 
group corresponding to an excited state of more than 400 MeV. excitation 
to about this fraction of the intensity of main group. The use ofa much 
higher bombarding voltage to increase the yield was impossible since the- 
alpha particles would have been obscured by protons scattered from the 


Fig. 4. 
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target backing of thick copper, and targets of the B isotope deposited 
on thin foils were not available. In order to extend the work to higher 
voltages a search was made for singly charged alpha particles corresponding 
to the excited state transition ; fig. 4 shows the results obtained with a 
target of thickness 0-165 mgm./cm.?, a proton energy of 430 keV., an 
angle of observation of 135°, and an analyser resolving power of 50. The 
analyser was calibrated at the higher energy values by observations of” 
the protons from the ?C(dp)8C reaction (Q=2-70 MeV.). The distribution 
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of particles observed with high analyser fields could not be completely 
resolved into separate groups, owing probably to the retardation of the 
recoil groups in the target. If, however, correction is made, as described 
in $2, for the effective depth of the target and for the presence of a 
carbon layer, the expected positions of the singly and doubly charged 
alpha particles and singly and doubly charged 7Be nuclei corresponding 
to the ground state transition are as marked (g, c, f, b). The groups (e, d), 
could be due to singly charged alpha particles and Be nuclei corresponding 
to a transition to an excited state of about 400 keV., and there is some 
slight evidence for an associated *Het++ group (a). In confirmation of 
this interpretation it was shown that all the groups except c and g 
(of energy 740 keV.) disappeared when the zinc sulphide screen was. 
covered by a 1-5 mm. foil, which indicated that groups a, b, d, e, f were 
of lower energy. Groups f and e were found to move with increasing: 
bombarding proton energy in the way expected of 7Be+ and 4Het nuclei 
respectively. 

At this time it was learned that definite evidence for the existence of 
a 435 keV. state in 7Be had been obtained by Brown, Chao, Fowler and 
Lauritsen* with better experimental conditions, and it was felt that 
nothing further could be gained by continuing work on the low energy 
group with the available targets. The evidence for the existence of an 
excited state of “Be near 400 keV. is in agreement with the results of 
Brown et al. ; the ratio of intensities of the transitions to the ground and 
excited states for a proton energy of 430 keV. cannot be estimated 
accurately from the present results, but appears to be about 5 to 1. 

The existence of a 435 keV. state in “Be should lead to the emission 
of two groups of neutrons in the reactions 


LAC) Ons a: meee tat oes. (3) 


and 
SUMO) Deas erin sabes ted an «1 (4) 


and to gamma radiation of energy 430 keV. in these two reactions and 
in the 1°B(pa)?Be reaction. Gamma radiation of this energy has been 
observed in the (pa) reaction and in reaction (4) by Lauritsen and 
Thomas *, and in the (px) reaction in this laboratory by Jelley and Freeman 
(unpublished). The observations of Freier, Lampi and Williams (1949). 
and of Taschek and Hemmendinger (1948) revealed no second neutron 
group in reaction (3) in intensity exceeding 10 per cent of the main group.. 
Recently, however, evidence for neutron groups in this reaction corre- 
sponding to levels in 7Be of energies 0-2060-07, 0-470+0-07 and 
0-745--0-07 has been obtained by Grosskreutz and Mather (1949). The 
745 keV. state would not have been observed in the present (pa): 
experiments, and the 470 keV. state seems to correspond to that found 


* We are indebted to Professor T. Lauritsen for sending us this information 
in advance of publication. 
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by Brown et al. and also indicated in the present work. Alpha particles 
corresponding to a transition to a 205 keV. state would have been 
observed in approximately the position occupied by group 6, which was 
interpreted as 7Bet++ recoils in fig. 4. The present work does not 
definitely establish the identity of group 6, but group f, which was 
interpreted as the corresponding *Be+ group, was shown to increase in 
energy with increasing bombarding energy as expected for particles of 
this mass and charge. 

Evidence for a second neutron group in reaction (4) has been obtained 
in the Cavendish Laboratory by Gibson (in course of publication) using the 
photographic plate technique. The neutron spectrum is consistent with an 
~ energy level of 435 keV. No evidence for a 745 keV. level was obtained 
by Mandeville, Swan and Snowdon (1949) in their study of this reaction. 


§ 3. COMPARISON OF THE REAcTIONS B(pa)‘Be, B(na)*Li. 


(a) Energy release. 

For the observation of reaction (2) the magnetic analyser and target 
-chamber were transported to the Atomic Energy Research Establishment, 
Harwell, and set up in such a way that a beam of thermal neutrons could 
pass through the 1B target. The region of the target of area 0-31 cm.* 
which had previously been bombarded by protons was defined by a slight 
discoloration and a screen was placed over the target for the (n«) 
experiments, so that only the alpha particles emitted from the area thus 
defined could be observed. A foil of stopping power 4 mm. was placed 
-over the zinc sulphide screen to prevent the detection of the 7Li+*+ nuclei 
produced in the reaction, while not appreciably reducing the efficiency 
for alpha-particle detection. The analyser resolving power was 25 and 
an angle of observation of 83° was used. The energy spectrum obtained 
is shown in fig. 5, curve (nx). The two alpha-particle groups, corresponding 
to transitions to the ground state and to the first excited state, are clearly 
‘evident. The main group in the (nx) experiment is appreciably broader | 
than that found in (pa) experiments, and this can be explained on the basis 
of target thickness effect by reference to a curve constructed by a procedure 
similar to that described for the derivation of fig. 2. In the (n«) case 
the cross-section for the reaction is the same at all depths in the target, 
so that the triangles representing the contribution to the total spectrum 
of the alpha-particles arising from various depths all have the same 
height, decreasing in width only with decreasing energy. The curve 
constructed on this basis is shown in fig. 6. It is seen that contributions 
to the peak are received from alpha particles emitted in the first 200 keV. 
equivalent thickness of the target, and the peak occurs at an energy 
110 keV. below that of alpha particles emitted from the surface. A 
‘comparison of this curve with that of fig. 5 shows that, in the position 
used, the target was effectively about 200 keV. ‘thick ’’, and that a 
-correction of 110 keV. should therefore be applied to the energy at the 
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peak of the spectrum. The energy of the main alpha-particle group 
emitted from the surface of the target was thus deduced ag 1-45--0-02 MeV. 
The effective thickness of the carbon film on the target was estimated, 
as previously described, to be 20 keV. for alpha particles, giving a true 
particle energy of 1-47+0-02 MeV. This gave the energy release in the 
MB(nax)’Li* reaction as Q=2-31+0-04 MeV. From the position of the 
alpha-particle group due to the ground state transition the true energy 
of these particles was estimated to be 1-78--0-03 MeV., giving an energy 
release of 2-80+-0-05 MeV. These values are in agreement with previous 
work (Boggild 1945, Gilbert 1948, Chao, Lauritsen and Tollestrup 1949). 
The accuracy of the present determination was limited by the intensity 
of the slow neutron beams available and by the solid angle of the analyser. 
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The 2°B target used in this (nx) experiment was bombarded with protons 
with the main object of comparing the yields from the two reactions, but 
the results also allowed another check on the energy release. A proton 
energy of 330 keV. was used, and an angle of observation of 83°. The 
energy spectrum (fig. 5, curve (pa)) indicated that the effective target 
thickness was greater than 120 keV., as expected from the (nx) result, 
so that a correction of 60 keV. for target thickness was necessary. An 
allowance of 20 keV. was made for the carbon film formed on the target, 
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and the true alpha-particle energy was deduced as 935+16 keV. This 
gave a value of 1-14-+0-03 MeV. for the energy release, in satisfactory 
agreement with the more accurate estimate given above. 


(b) Cross-sections. 

The ordinates of fig. 5 give directly the yield of alpha particles in 
the reactions (1) and (2) per incident proton or thermal neutron respectively. 
The yield for the proton reaction was obtained by integrating the charge 
collected on the target during the (px) experiments, care being taken 
to eliminate the effects of secondary electrons. The amplifier gain and 
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bias setting were known (see I., § 2 (c)) to be suitable for the detection of 
alpha particles of 800 keV. with 100 per cent efficiency, and were high 
enough to make correction for the presence of the lower energy 7Be++ 
nuclei unnecessary. In the (nx) experiment the flux of neutrons 
incident on the target was measured by placing a calibrated manganese 
foil of area 1 cm.* in the position of the target and allowing it to be 
irradiated for a known period. The beta-activityy was found in a standard 
counting arrangement, whose efficiency was checked by means of a uranium 
source ; the neutron flux used was 6-5 x 10°+5 per cent. No correction 
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for the presence of 7Li recoils in reaction (2) was necessary, since an 
absorbing foil was used over the zinc sulphide screen in the (na) 
experiment. 

The analyser resolution in both (nx) and (px) experiments was less 
than the equivalent thickness of the target. The ratio of yields must 
therefore be corrected by a factor (1-57) equal to the ratio of the energies 
of the particle groups. <A further correction for charge exchange of 
1-27 for the (nx) and 1-66 for the (pa) reaction was deduced from the 
figures given by Briggs (1927). With these corrections, and after 
subtraction of background counts, the ratio of yields at 83° per incident 
particle becomes 1-37x10-°. This, assuming isotropic angular distributions, 
is taken to be the ratio of the cross-section for reaction (1) with 330 keV. 


Bige7, 
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Reaction 1°B(n«)?Be excitation function. Absolute cross-section at 
330 keV .=<4:5 x 10-27 cm.? 
protons to that of reaction (2) with thermal neutrons. The neutron 
flux measurement by means of the 1/v-law manganese detector assumed 
a neutron velocity of 25x 10° cm. per sec., and the B cross-section for 
neutrons of this velocity is 3285 barns, from the results of Fermi, Marshall 
and Marshall (1947). The (pa) cross-section is, therefore, 
Oy (3380 keV.)=4-5 x 10-?” cm.? 

In order to compare this figure with that previously given for a proton 
energy of 530 keV. in I., an extension of the excitation function given 
in I., fig. 11, was made. This is given in fig. 7, the ordinate scale being 
fixed by the 330 keV. determination, and from this 
py (530 keV.)=1-9 X 10-?8 cm.?, 

Zee, 


348 Emission of Short-range Alpha Particles under Proton Bombardment 


which is to be compared with the value of 1-2 x 10-6 cm.? giveninI. The 
difference is not surprising in view of the uncertainty of the factors which 
entered into the earlier value ; in the present method knowledge of the 
solid angle of the analyser is not required, and the exact target thickness. 
need not be known since the analyser selects a definite layer of the target 
determined by the resolution. The calculation of the cross section given 
in I. from the yields shown as ordinates in I., fig. 11, included a large 
correction for the ratio of areas of the proton beam and the boron target; 
this correction is unnecessary in the present experiment. The main 
source of experimental error, which has not been eliminated in the present 
cross-section determination, is a possible over-estimate of the counting 
efficiency for alpha particles of 3-5 mm. range emerging from the stopping 
foil in the (nx) reaction. From the curves given in L., fig. 4 (b) (bias 20), 
this error is unlikely to exceed 20 per cent. 

Asa further check the cross-section for the (px) reaction was also deduced 
from the observations on small weighed targets described in §2. The 
solid angle of the analyser and the target depth determined by:resolution 
were assumed and a value of 2-0 x 10-76 cm.? at KE, =530 keV. was found, 
in agreement with that deduced from the (nx) reaction within the limits. 
of error of the measurements. 
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SUMMARY. 


Orthodox analysis of most problems of plastic flow is extremely difficult, 
in that the governing equations are not linear, so that partial solutions 
may not be superposed. Many problems, however, yield to treatment 
by relaxation methods. This paper considers the consequence of a 
particular hypothesis regarding plastic strain, viz. the Mises-Hencky 
eriterion, and, by application of a numerical method due to Allen and 
‘Southwell, traces the effect of plastic straining. The success of the 
method in no way depends upon the particular yield hypothesis adopted. 
‘Three problems are considered: 

(i) A two-dimensional slit tension member. 

(ii) A sharp-ended bar under specified edge tractions. 

(iii) The uniform compression of a block between two smooth, rigid 

plates which do not overlap it. 

In each case, results and conclusions are summarized and presented 


‘graphically. 


§ 1. INTRODUCTION. 


ORTHODOX analysis of most problems of plastic flow is extremely difficult 
in that the governing equations are not linear, so that partial solutions 
may not be superposed. Moreover, in many cases the shape of the 
boundary and the boundary conditions are such as to preclude an exact 
solution, even in the elastic case. Many of these difficulties are resolved 
when the problem is attacked by relaxation methods—thus, for example, 
the shape of the boundary does not materially increase the complexity 
of the problem. Again, many of the exact solutions are extremely 
complicated, being usually expressed in series form. The convergence 
is often very slow, and the numerical value of the stresses at various 
points would require as much labour as a normal solution by relaxation 
amethods. 

Relaxation methods are peculiarly suited to a class of problem involving 
boundaries or interfaces not initially known—problems which, in many 
‘cases, present insuperable difficulties to an orthodox approach. In this 


* This paper is based upon work forming part of the writer’s Ph.D. thesis, 
‘London University, June 1949. ; 
+ Communicated by the Author. 
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connection, besides problems of plastic flow, relaxation methods have 
been successfully applied by Shaw and Southwell (1941) to problems 
relating to the percolation of fluids through porous materials, and by 
Southwell and Vaisey (1946) to fluid motions characterized by free 
streamlines. 

The true criterion of elastic failure and the relation between stress. 
and strain during plastic distortion in real materials has yet to be 
determined. Many hypotheses have been advanced, the validity of 
any of which only experiment can finally endorse. Direct experimental 
study of stress-strain relations is made more difficult by the fact that 
stress is not itself measurable. Experiment can only verify or deny 
results computed from one particular theory in regard to total displacg- 
ments produced by imposed external tractions. We adopt a particular 
hypothesis regarding plastic failure and use it as a basis for our 
computations. It is to be emphasized that the methods evolved do not 
depend on the particular hypothesis adopted, but could be applied 
similarly with any other criterion. 


Fig. 1. 
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§2. GENERAL THEORY. 


The first attack on plastic problems by relaxation methods was made. 
by Christopherson and Southwell (1938), who considered the plastic 
torsion of a bar of equilateral triangular cross-section. Later 
Christopherson (1940) applied the same methods to I-section girders, 
and Shaw (1942) considered hollow (i.e. multiply-connected) sections, 
where plastic strain can commence at the external or at an internal 
boundary. All these problems of plastic torsion based their solution 
on a hypothesis put forward by Prandtl (1923) which assumes the relation 
between shear stress and shear strain to be that indicated in say ed Be 

For stresses below some limiting value f, the material behaves 
elastically, but the stress cannot exceed this limiting value f, with which 
any shear strain in excess of a limiting elastic shear strain Yy can be 
associated. This imposes (in the torsion problem) an upper limit upon. 
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the gradient of the stress function. This restriction can be represented 
by what is usually known as the Prandtl roof, and the ensuing computation 
is straightforward. 

The next major advance was made by Allen and Southwell (1950), 
who investigated more complex problems of plastic straining in two- 
dimensional stress systems. In particular, they considered notched 
tensile specimens having a central waist. Two examples were taken, 
in which the notches were semicircular and sharp re-entrant right angles. 
One of the examples we have considered is a limiting case of this Ey pe 
of problem—a slit tension member. 

For the more complex stress systems considered, a generalization of 
the Prandtl hypothesis is required. The generalization we adopt rests 
on two assumptions : 

(i) The first occurrence of plastic straining is determined by the 
** Mises-Hencky ” criterion, 7.e. the normal elastic relations hold until 
the principal stresses p,, po, p, attain such values that 

(P2—Ds)*+(P3—P1)*+(Pi- Po? = 6h, « « . « (I) 
k being a constant for any material. The relation (1) holds thereafter 
throughout any region of plastic strain. 

(ii) Where the material has yielded plastically, the strain consists 
of both an elastic and a plastic component. For the elastic part the 
usual stress-strain relations hold, viz., 


Ci 5 {P:—2 O(PoaieDa)ime-= eee! vo fe, + 22) 


etc., where E (Young’s modulus) and o (Poisson’s ratio) have normal 
values. 
The plastic part is assumed to be equivoluminal, 7. e. 

ep eye, a0 ee renee Se a. (38) 
and in it the incremental shear strains resulting from the three principal 
shear stresses are directly proportional to these stresses and occur on 
the same three planes, 7. ¢. if 71, ¢2, ¢3 denote the prinicpal shear stresses, 
and Ay,,Ay2, Ay; the principal incremental shear strains, then 


Ay: AY: : AY3=%1 : V2: Is- SMa iees. tac (4) 
Then if Ae,, Aes, Ae, are the principal incremental plastic strains, 
(4) may be written 
(Ae,—Ae;) : (Ae;— Ae,):(Ae,—Aes)=(p2—Ps) :(Ps—P1) ?(Pi—P2)- (5) 
where Ae,, Ae, Ae; have the same directions as p;, P2, P3- Equations 
(4) and (5) may be written 

eS en een oe ny, .',".'(6) 

q1 q2 qs 


and 
(Ae,—Ae,) _ (Ae;—Ae,)_ (Ae; ACs) _ gy) > (7) 


Do—Ps (is PP Py—P» 


bol 
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where Ad is an incremental quantity. 
‘Thus 2Ae,—Ae,—A0e,=3 4X(2p,—po—Ps), 
i. é€., using (3), 
Ae,=AA{p,—#(potps)}, - »- - - + + (8) 


Thus, as regards the plastic component of the distortion, the stress- 
incremental strain relations have normal sees the effective values of 


Young’s modulus and Poisson’s ratio being —. and $ respectively. 


a 
The total strains associated with given stresses are given by 
“€y=e, 7 ey 


etc. That an exact analysis is rendered difficult is shown from the fact 
that (2) relates to total and (8) to incremental strains, both expressible 
in terms of total stresses. 

Combining (2) and (8) gives three equations of the type 


Ae,= = {4p,— o(Ap.+ps)}+ 4A{pi— 3(Ds-+-P3)j-0 ed 


Allen and Southwell (1950) deduce from the Mises-Hencky criterion (1) 
that for plane stress 


“Xie XY V24- 3X2 <3, fog ee ee 
whilst for plane strain 
(X,—Y,)2+4k2<4k, 2 2 2. 2 2 (il) 


the sign of inequality relating to elastic conditions and the sign of 
equality to conditions of plastic straining. The condition (10) follows 
immediately, but (11) is obtained on the simplifying assumption that 


o=%, whatever value may be appropriate in (9). The error due to this 


assumption is small, 


§3. THE COMPUTATION OF THE STRESSES. 


In plane strain (and in plane stress), provided that there are no body 


forces, all of the conditions of stress-equilibrium may be satisfied by 
writing 


7 Oy? ~Y Oxdy? ~Y Oat ny 


whatever be the nature of the accompanying strains. When known 


tractions act at a single boundary, y and have the same boundary 
Vv 


values whether plastic straining has occurred or not. 
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For specified edge tractions N,, N, we have that 
N, cos (x, v)—N, sin (x, v) =X, 
=X, cos (x, v)+ X, sin (2, 


v) 
ty 0 d)\0x a (ox 
= {cos (xv, v ay — sin (a, ates a aa F, (3%). 


In a similar manner 
N, sin (x, v)+N, cos (a, v)=Y, 


=Y,, sin (x, v)+ X, cos (x, v)=— g (2). 
Hence for all points of the boundary 
os ox : 
=[¥, de, HA = [Xa 22) eeeiarees (13) 
both quantities being single-valued. 


Fig. 2. 


: vy 
<2 y Ns Ny 
; Xy 
ac 


The specified boundary tractions can thus be reduced to given boundary 


5% 


0 : 
values of a sk. Since the stress components depend on the 
-derivatives of ey st, the lower limits of integration in (13) are 
: Ox dx 
arbitrary (so that Dae may be given the value zero at some 


-convenient point of the boundary). 
Finally y can be calculated at all points of the boundary from the 


‘equation 
Ox dx 4 Ox dy 
x= {ee 7a au eh as, ru ey oe (Le) 


the lower limit again being arbitrary. 
Writing (14) in the form 


Loy ae ey 0x ox 

x" Ox gue "ash 5k hes 

a oe He Keath? Werte kb 
aa v3 {@x,— —yX,) ds, (15) 


it follows that y is single-valued at all points of the boundary. 
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) 
Equations (13) and (15) enable boundary values of fe aa to be 


computed from specified edge tractions. So long as the whole of the 
plate remains elastic, the determination of x at internal points is a 
biharmonic problem, whose method of solution by relaxation methods 
is given by Fox and Southwell (1945). Apart from a multiplying factor 
which increases linearly with the total load, the solution is unique. By 
substitution in the criterion (11) a value can be found for the greatest 
load that can be sustained elastically. If the load is increased still 
further, the criterion will be violated in some regions. In such regions 
the elastic solution must be modified to satisfy it, whilst leaving the values 


0 
of x, =, 3 unchanged on the boundary. This modification~ can 
be carried out by a relaxation method, described by Allen and Southwell 
(1950). 


eeece 1 2 os SS ALS > x 


§4. PROBLEM (i.) SrRESS DIstRIBUTION IN A TWO-DIMENSIONAL 
Suir TENSION SPECIMEN. 


Fig. 3 shows the slit tension specimen, subject to specified edge tensions. 
Three cases are considered, in which the ratio qT of the total width of 
the bar to the width at the waist is given the values 8, 4 and 2. From 
symmetry it will be sufficient to consider the quadrant OABCD. We can 
eliminate dimensional factors by writing 


ta Dy Y=Dy, Vee Ly! heges Wael oa 


2T denoting the total tension on the specimen. 


/ 


saa x ‘ ft 
Then such quantities as =4 measure stresses as multiples of =, 


Oa’? 
the average tensile stress at the ends. Hereafter, however, the dashes. 


attached to x, y and y will be suppressed, and the symbols will stand for 
purely numerical quantities. 
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The boundary conditions are, along AB, 


at 


nex 
| ‘ (17) 
is specified (and was taken as 3200 to obviate decimals). 
ewe —y - 
Gedy” 
along BC and CD, X,=Y,=0, i.e. 
0 (ox 0 (dx 
—(——)=~(~—)= i i-cth Leese LS) 
(3) ds (Fr) : Sod 
so that ox and ox are both constant. 
Ox oy 
If we take He =0 at D, then x _9 along DC and CB, and hence 
ox oeecae B d is th lt 
Oa =0 along AB (since az 8 constant along BA and is thus equal to 


its value at B). 


: 07x Ox tO Xas ral Saas 
Since cn eee along AB, Bs =3200y, taking ay =0 at A. us 
at B, (y=1), St — 3200, and hence a =3200 along BC and CD. 
Further, taking y=0 at A, 
N= 1600; "04 = ae eens car (19), 
along AB. Finally, from (14), 
Ox dx ne OX a} as 
See i Ox" as oy ds 
_ (ex dy 
— Agadaee 
along BCD where ox 2-0 
Ox 
= | 3200 dy 
== S200Ys= GOD ee Pay oie eS als (20); 


since y=1600 at B (y=1). Actually along BC (y=1), x=1600, and along 
CD, x=3200y—1600. 


The case ?=4 was first attacked. After a preliminary investigation 


on a coarse mesh, a finer net (h=4) was used. 
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The computation of y and the stresses in the elastic case was carried 
out by the relaxational methods given by Fox and Southwell (1945). 
From this elastic solution we have to determine what modifications are 
necessary when, in a region of plastic straining, the condition (11) super- 
venes. It appears that a plastic enclave * will start at the bottom of the 
slit and extend into the specimen ; within an enclave, the initial (elastic) 
distribution of y has to be modified in accordance with (11). Allen and 
Southwell (1950) give a relaxational method for the computation of x 
within a plastic enclave. The process is rather laborious in that the 
relaxation pattern differs from point to point, and strictly only applies 
to infinitesimal increments in y. Much labour can be saved by the use 
of inexact patterns, exact account, of course, being kept of the residuals, 


eal 
Y 


=) 


Values of the Mises-Hencky equivalent stress (taking 1000 as its maximum 


value) when the plastic enclave has spread right across the bar. Cases =4. 


In all the examples considered here, a modified relaxation technique was 
adopted—inexact patterns being used to give the general trend of the 
x distribution, the final values being obtained by a systematic trial-and- 
error method. The labour involved was not prohibitive. The determina- 
tion of the growth of the plastic enclave is facilitated by the computational 
device adopted of proceeding from step to step by lowering the value of 
the maximum equivalent stress rather than by increasing the applied 
tension. 
ne eee en MN Se ea LA el 

* Southwell has suggested the use of the word enclave to denote a region 
of plastic straining. 
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Four distinct stages in the growth of the plastic region were obtained. 
Fig. 4 records the values of the Mises-Hencky equivalent stress when the 
plastic enclave has spread until it extends right across the bar. Fig. 5 
summarizes these four stages, and shows the growth of the enclave with 
increasing tension. 


, the ratio of 


The same problem was considered for other values of a 
a 


D 
the total width of the bar to the width at the waist. The case a 7 


ne a D 
exhibited the same general characteristics as the case 7 =4 already 


Fig. 5. 


Diagram showing the progress of plastic straining—superposed plastic-elastic 


D 
interfaces. Case — =4. 


d 


considered. The elastic solution was verified by photo-elastic methods, 
and excellent agreement was obtained in spite of the difficulty of obtaining 
a model with a “slit ’—practical considerations necessitating the notch 

: D 
having a finite though small width. One difference in the case 7q 252 
is the rapidity with which the plastic enclave spreads. Although the 
growth of the enclave exhibits the same characteristics as in the case 


4, only two stages in the growth were obtained, as compared with 


Cee, 


; . . q F 
four in the previous case. Once plastic flow begins, the enclave spreads 
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-extremely rapidly and extends right across the bar for a very small 
increase intension. Fig. 6 gives the values of the Mises-Hencky equivalent 
stress when the plastic enclave has spread right across the bar. 


Fig. 6. 


Values of the Mises-Hencky equivalent stress (taking 1000 as its maximum 


value) when the plastic enclave has spread right across the bar. Casey =2. 


D 
The case ae =8 was an attempt to solve the problem for an infinite 


plate. There was a marked distinction in this case as compared with the 
other two cases. As plastic flow begins, an enclave, of the same shape as 
before, forms at the root of each notch and spreads into the plate ; this 


Fig. 7. 


is accompanied a little later by two new enclaves starting from points 
H, H on the centre line near the waist (see fig. 7). This phenomenon 
was also found by Allen and Southwell (1950) in the case of notches 
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which were semicircular. The enclaves extend extremely rapidly and, 
: D 
as in the case 7d =2, only one intermediate stage was obtained before 


the two enclaves coalesced. Fig. 8 gives the value of the Mises-Hencky 
equivalent stress for the most advanced stage of plastic straining obtained. 


ae) 
As the ratio — increases, so the labour of computation greatly increases, 


d 
because the number of mesh points rapidly multiplies. The case 5 =8 
Fig. 8. 
A a a A a zsol|_~ 29] anal 368! 330] __sss|_ =| 334 325] 3344 _ 
57 ya} tay zag] sol _ssit_sai| ai aa ar | sg] -35c} 
B31 219 seal ser} seal 438 4511 460) 438 ais] 38a) 385| 382) 
12] 292] 438] 47a} sul 4uol 42a] __ 398 22] 
[202] __saa| s3zl 553) 5781 590] Syl 5351 a 4021 Ha! 
A Lie, eee agen S| tse 0 ssl seal fst aul 
34g 65 Gaol ges} so|___esa_ est sel __s aca 27 
oS} __7é0} can G51) is} 355} Seal 
PA 
i i mall 
| a} 
ae ---aall_Lohd \__ eas] 


& 


Values of the Mises-Hencky equivalent stress (taking 1000 as its maximum 


D 
value) for the final plastic stage. Case 7q =o. 


should give a good indication of the results for an infinite plane, as there 
was little variation in the stress distribution a long way from the slit. 


D : a 
Preliminary investigation for the case = =10 showed but slight variation 


d 


from the case = =8, 


§5. SumMARyY oF RESULTS. 
Hill (1949) considers the state of stress in the core of a notched bar 
under an applied load sufficient to cause plastic flow, the notch being 
of sufficient depth to ensure that the plastic region is localized in and 
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around the core. If the curvature of the notch is greatest at the root,. 
the yield limit is first attained at that point, and as the applied end load. 
is further increased, Hill assumes that the plastic region will spread in 
some manner from the notch to the axis, and shows that this plastic 
region must be of a certain minimum extent when the bar first extends. 
The initial growth of the plastic region for the limiting case of a deep 
notch, viz. a slit, which has been considered here, is in a direction which 
could hardly have been forecast in advance, and is not towards the axis 
of the bar as suggested by Hill’s diagrammatic figure. Nothing in 
the present paper, however, conflicts with his results and conclusions :-— 
“So long as part of the core is not plastic, no displacement of the 
ends of the specimen is possible. Even when the plastic region has 
spread right across the core, it does not necessarily follow that deformation 
is possible, since the plastic material may be rigidly constrained by the 
adjacent bulk of non-plastic material’. From a consideration of the 
equations of Saint Venant governing the flow velocities in the plastic 
region, Hill shows that no extension of the specimen is possible until 
the plastic region has spread sufficiently far to include the entire slip 
lines OS, OS’ from the centre O to the free surface of the notch. At the 
moment when extension first occurs, the plastic boundary must touch 
OS, OS’, or must pass through either S and 8’, or through O. Hill’s 
surmise that the plastic boundary actually passes through 8 and §’ is 
confirmed by the results of this work :— 

‘““ Many writers, including some of the most recent, such as Sokolowsky 
(1946), and van Iterson (1947), have suggested slip-line fields for various 
problems quite arbitrarily, with no adequate discussion of their validity. 
In some cases elementary considerations are sufficient to prove them 
wrong, but in others only detailed numerical calculations can decide. 
In many problems the slip-line field in the region of plastically deforming 
material may be surmised without a full solution for the stresses in 
both elastic and plastic regions”’. The work of Allen and Southwell 
(1950) and the results of the investigations presented here show that 
in many cases the actual plastic region develops in a direction which 
could hardly have been surmised in advance. Thus, although it may 
be easy to suggest apparently plausible slip-line fields for a number of 
problems, detailed numerical calculations constitute the only certain 
method of procedure. In a future paper the results will be given of the 
investigation of a sharp-ended bar under specified edge tractions and of 
the uniform compression of a block between two smooth, rigid plates. 
which do not overlap it. 


§6. ADDENDUM. 


Since the above investigations were carried out, Hill has shown (1950) 
that in the plane deformation of an incompressible plastic-elastic material 
under a given loading-path the distribution of stress at any stage can 
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only be found uniquely from the stress equations alone if the plastic- 
elastic boundary is not cut more than once by any slip-line (and if, 
further, unloading does not occur). His proof depends on the fact that 
the slip-lines in the plastic region are the characteristics of both the 
stress-rate equations and the velocity equations. It follows that 
the calculated plastic regions in this paper are, after a certain stage, 
possibly unreliable since theoretically an unlimited number could be found 
to satisfy the stress equations and boundary conditions. Uniqueness 
is obtained by choosing the plastic region so that a velocity solution 
can be found. The results given here, however, do exhibit correctly 
the growth of the plastic region in the initial stages, and it is not 
anticipated that the final results are much in error. 
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XXIX. On Tolman’s Equation describing the Thermal Equilibrium in a 
Gravitating Sphere of Fluid. 


By H. A. BucHDAHL, 
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[Received November 9, 1949.] 


ABSTRACT. 

It is shown that by adopting a form of the line element slightly 
different from that chosen by Tolman a first order equation for the 
dependence of proper temperature on gravitational potential in a 
gravitating fluid sphere appears in place of Tolman’s second order 
equation. The problem suggested by him concerning the meaning of 
a certain constant of integration therefore does not arise. 


. $1. INTRODUCTION. 
Tortman (1934) has investigated the distribution of temperature in a 
gravitating sphere of fluid by invoking the thermodynamic conditions 
for static equilibrium. Assuming the form of the line element to be 


ds*= — exp {u(r)}(dr?-+-7r? d0?-+-7? sin? 6 dd?)+ exp {r(r)} d?,. (1.1) 
he arrives at a differential equation of the second order, a first integral 
of which is 

dlogT,  _ldv_ Bexp{—}u—}}T, As 
gp Bg he ee 


where B is a constant of integration and Ty is the absolute temperature 
as measured by a local observer at rest. In order to arrive at a solution 
of (1.2) in the form 


T, exp {iv} =T,V/(gaa) = const. “51 eed tee 


he sets B equal to zero “on physical grounds” and remarks that “ it 
would of course be interesting also to investigate the possibility for 
solutions of physical interest in which the constant of integration was 
not taken to be zero ”’. 

We shall show, however, that the appearance of the term involving 
B in (1.2) is spurious in the sense that by adopting a form of the line 
element slightly different from (1.1) there arises a first order equation 
which replaces Tolman’s second order equation, and which has the integral 
(1.3). The question of having to assign a particular value to a constant 
of integration equivalent to B consequently never arises. It might be 
added that we are here not concerned with the question of the validity 
of Tolman’s thermodynamics; and accordingly we make full use of his 
method of attacking the problem in hand. 


* Communicated by the Author. 
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§2. DERIVATION OF THE EQuation. 


We adopt a spherically symmetric line element in the form 
ds*==— exp {X(r)} dr?—r?(d6?+-sin?@ d¢2)-+ exp {v(r)} dt. . (2.1) 


In the following a number preceded by the letter T refers to the equations 
of the work of Tolman quoted above; and his notation is adhered to 
throughout. Then, by T(95.13), 


Ht GT! 1 
8apy= exp (-a(= +3) Sieg: itch nen? i g(22) 
3 fA 1 


‘Consequently 

5H =8(4arpo9 exp {4A}? dr) 

== (—2rpogr? +4) exp {4A}SA dr+4r exp {—4A} 80’ dr, 
by:(2.3). Also, by (2.3), 

Pop =P o9(47 exp {FA}? dr) = 2rpor* exp {FASS dr. 
Hence 
OB) + 199Vp=[—277( Pon — Po)? +4] exp {ZA}OA dr+4r exp {—4$A}8)’ dr 
z = {[4r(v’—A’)+d]dA+4rdr'fexp{—tA} dr. . . . . (2.4) 

Inserting (2.4) in the condition for thermodynamic equilibrium T(124.5) 


we have 


pees 
ee (yro—N) FAFA} dr—0 . . . (2.5) 


for all arbitrary variations which vanish on the boundary of the region 
considered. Integrating the term containing 6A’ by parts, and noting 
that the integrated part vanishes by hypothesis, (2.5) becomes 


ESD ios} oy ey rexp{—}\’ 
{Se ae Fe gees (“") }oa dr 


, l , 
ee ay abo OLR a 4 
=f exp { | sar @ Joa dr=0. 


Accordingly we must have 
TN Ge ; 
2 Nee) 
(7,) —3n, sy 


_'This is the first order equation which replaces T(128.2); and its integral 
is again given by (1.3). 
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XXX. Scattering of Pseudoscalar Charged Mesons by Nucleons.—IL. 
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SUMMARY. 
In §1 the fourth- and sixth-order contributions to the scattering 
cross-section are evaluated for the non-relativistic region (w+3y). The 
other section deals with damping effects. 


§ 1. SIXTH-ORDER CORRECTIONS TO THE SCATTERING CROSS-SECTION. 
WE shall now investigate the same fourth-order transitions which were 
discussed in §5 of Part I. (Corinaldesi and Field 1949). A more refined, 
still non-relativistic, treatment will be slat ae for 

Pi PE PF he <i (1) 
- ww Bw ial f Ay) 5 
but with 
Pe Pr iy Pr™e, (w > 3p). . . . . . . (2) 

In order to simplify the calculation we shall choose as our reference 
system that in which the centre of gravity of the proton and meson is 
at rest. This will not affect the value which will be obtained for the 
cross-section. In the new reference system, 


wal wei ane ee 


by ae yD 
From the conservation of energy between the initial and final states, 
it follows that 
[Pil =[ pel =[Pyl =Lp;| (=p, say) eG) 
and 
(Py « P7)=p* 008'0,— 0.0 eve oe ee 
where @ is the angle of scattering. We have also, with the notation 
of Part I., 
W;=W p= 
and 


2= 2,=2. . . . . . . . . (6) 


* Now at Physics Department, University of Manchester. 
+ Now at Department of Mathematical Phy sics, University of Birmingham. 
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When conditions (1) are (2) are imposed, only eleven transitions are 
found to give contributions of the fourth-order which are not negligibly 
small. They are the same transitions that were discussed in Part I. 
All others have to be neglected, even for the higher energy range now 
under consideration. 

Analogously to formula (49) of Part I. we have now 


Wot, >= ap 8 { {(#,(Q+a.P,+MB)H# ,(Q+a.P,+MB) 
+ # (2-+-a.P:+Mp)#,(Q+-a.P,+Ms)} . . (7) 
where 
_ 2g? fw—Q+0. LMA ; 
ai To | (Q—woP— OJ" 8 
H OS ee er {Pys{ 2(P3) 
* 2N?w y o(p’).2(p;) 2(p2) 2(p;)(E;,—E)) (E;— E,,,)(E;—E,,) : 
Fa. p,tMB)By,| 2(p.)+o . po FMP ]By5[ 2(ps)F-@ - pst MB IBy5}- 
(9) 
In (9) 
P2=P’. Hh 
pPi=p +p; or —P, Le ee Re 1O) 


Ps=p'+p; or —P. 
Performing the integration over momentum space, and neglecting terms 
of second and higher orders in p/M as small compared with 1, we obtain 


37 Sts 8 ao) CD eer ul 
Ha, HH, 92 pate (243g 3 oO) . ° ( ) 


Similarly, | H, |? is evaluated from (8) and formula (16) of Part I. 


pes mgt op S ten eles 
P= ype (1+ i eo): oe (12) 


The differential cross-section, do, is then given by 


do=do,+do,+ . 


eee, ay eg .~<|,. + -» <<,(13) 
(2)? 
in analogy with formula (55) of Part I. Thus the fourth-order term is 
dow= Tap ae ee i a) sin 8.8 Ole ee. Beet (4) 


and the sixth-order term 
‘ Rone gee 0 dod ait 
deg—=— leche a (2+35 = 8 a) sin 6 dé dd, (15) 


where sin 6 dé d¢ is the element of solid angle. 
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Integrating over @ and 4, we obtain for the fourth-order contribution: 
to the total cross-section 


7 w 
Cpie= me (+s) a dae ee, woe (14 ) 
and for the sixth-order contribution 
3 w 
%%=— Te” (2+33). « oh rs ° (15 ) 
The ratio of the sixth- to the fourth-order terms is. then given by 
ae Sok a 16) 
alae (2+3)- - “yee (16), 


This formula shows that, to the order of approximation used here, the 
ratio of the two contributions to the cross-section increases linearly with 
the energy. As for the convergence of the series o,+o5+ ..:, this. 
depends, as was pointed out in Part I., on the choice of the numerical 
value of the coupling constant g. 

At the present stage of the development of meson theory it is too- 
early to say which of the values proposed by various authors is likely 
to be in agreement with experimental results—in so far as these already 
exist. In fact it is doubtful whether a pure pseudoscalar theory with 
a pseudoscalar coupling is likely to give a satisfactory account of the 
scattering process considered here. Several authors have adopted the 
value g?~36, which would make the series divergent (Villars 1947, 
Dyson 1948, Luttinger 1948). Slotnick and Heitler (1949) have suggested. 
a slightly lower value, g2~20. On the other hand, Moller and Rosenfeld 
(1940), in their mixed theory, have chosen for the pseudoscalar coupling 
of pseudoscalar mesons a value g?~0-1, for which the series would 
converge. Recently, Bethe (1949) has been proposing yet another value, 
g~4., 

Finally, we want to point out that the numerical value obtained for og 
depends on the way in which the renormalization of g is carried out. 
Since this point was not discussed in detail in Part I. of this paper, we 
shall give a short summary of the procedure here* ; 

After having carried out the mass and “charge ”’ renormalization, 
we aré left with og=og+-0%. Here, og is what remains of the previously 
divergent matrix elements, and og is convergent. Now, og vanishes in 
the extreme non-relativistic case (E.N.R.) and even in the approximation (1) 
considered in this paper it does not give any contribution to the cross- 
section (15’), which is due only to og. In the case of the Compton effect, 
as treated by Corinaldesi and Jost (1948), the situation was different. 
There og, too, was zero in the E.N.R. case and this circumstance made. 
it possible to say that the renormalization of the electric charge had 


Sea TOT neh? ene s Dw k teh o 
* We are grateful to Dr. M. P. Schafroth, Ziirich, for asking us to clarify 
this point. 7 
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been performed in such a way that limo,=0 for p>0. If we would 
apply the same criterion in the present calculation, the whole contribution 
to og which is constant for p+0 would have to be left out, and og would 
be much smaller than o, for all energies under consideration. But we 
see no reason why we should apply this criterion to meson processes, 


§ 2. Dampine EFrFrEcts. 


Another correction to the cross-section may be obtained if we consider 
the reaction of the emitted meson on the motion of the nucleon. Attempts 
have been made by several authors to take account of this damping 
reaction (Heitler and Peng 1942, Wilson 1941, Sokolov 1941, Gora 1943) 
and we shall briefly discuss its effects on the results obtained above. 

The scattering cross-section may be written in the form 


op=27p(w)| Ul, Sei oli gic eee Aes scl) 
where U is a probability amplitude defined by the integral equation 
U=K—i7KpU, ree ese “eter 8) 


U is a matrix element depending on the initial and final states, K is a 
similar matrix element given by K=H,+H,+ ... and KpU is a matrix 
product, with summation over all final states. In order to obtain an 
idea of the order of magnitude of the damping effect it will not be necessary 
for us to solve the integral equation (18) rigorously. Substitution for U, 
in the damping term, of an average over angles will be sufficient for our 
purpose. We then obtain 


K 
DS 19 ] 
he 1+impK’ 2) 
so that the corrected cross-section op is given by 
[ey 
2 ee 20 
oH 1+ 2) ( ) 
where 
j Ge) | eee se (21) 
and. 
fon p*| 1 |? Nee er ee te clan) 


We shall now discuss the behaviour of oy for different values of g?. For 
very high values of the coupling constant the series for |KP diverges. 
mp 

but there is hardly any sense in including higher-order terms in the 
expressions for |K[? and «?, and this result would be very unreliable. 
In this case, we might rather assume that the first non-vanishing approxi- 
mation of K comes closer to describing the real physical process than the 
inclusion of the rapidly increasing higher-order approximations would. 
Then we should have to put |K|?=|H,[?. At very low energies, the 


The damping would still make the cross-section finite, namely op=— , 
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values of op for g?=20 are, of course, very much higher than those for 
g?=4. With increasing energy, however, damping is much greater in 
the former case (since «2 cg‘) and at the upper limit of the range for which 
the formule are still valid (see (2)), cross-sections for both values of g* 
are of the same order of magnitude (~2 x 10~*6 cm.®, see fig. 1). 


Fig. 1. 


25 


is) 
° 


7 . 
Ig*x10"Cm3—> 


th Ooh 


2u ym 4u “ 6% 


A 


Energy dependence of damped cross-section, with | K ?=|H,|?. 
(i.) o/g* (no damping), 
(ii.), (iii.) and (iv.) op/g* for g?=1-2, 4 and 20 respectively, 


For lower values of g*, ¢. g. g?=1-2, the damping effect is much smaller 
than in the previous cases. Moreover, it is now reasonable to include the 
sixth-order terms in the expression for |K[?. Fig. 2 shows the damping 

7 : Ci be 2. pp 2 j 
effect for the two cases | K =| H,|? and| KP=|H,.?+H.4+H, ». 


In conclusion, we would like to thank Professor W. Heitler for much 
valuable criticism. 


Scattering of Pseudoscalar Charged Mesons by Nucleons 
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Influence of sixth-order correction on cross-section, for g?=1-2, 
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(iv.) op/g* with |KP=|H, P+H, a+ Hg». 
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SUMMARY. 


Some recent studies have indicated that in molecular-orbital calculations. 
on alternant hydrocarbons the explicit inclusion of the overlap integral, 8, 
between adjacent atomic orbitals makes no significant difference in the 
consistency of: the results (Wheland 1942, Brooks 1941, Coulson and 
Chirgwin, in press). In the present paper a further investigation is 
reported on total z-electronic energies and the various polarizabilities of 
such hydrocarbons. In connection with a regularity observed by 
Wheland (1941), the resonance energies are discussed in relation to the 
total z-electronic energies ; it appears that the latter are more consistent 
than the former. A correlation between the atom—atom—‘ self-polariz- 
abilities’ and the conjugation energies across single bonds leads us to 
define ‘‘ conjugating powers” of hydrocarbon residues in complete 
concordance with the analogous definitions given in the simpler theory, 
neglecting overlap. 


INTRODUCTION. 


ONE of. the most surprising features of the molecular-orbital method, if 
applied to alternant hydrocarbons (for a definition, see Coulson and 
Longuet-Higgins 1947 b), is the fact that the explicit inclusion of the 
overlap integral, 8, between adjacent carbon 2p.-atomic orbitals makes 
no significant difference in the self-consistency of the calculations of 
several molecular properties. This was first shown by Wheland (1941) 
for the resonance energies, by comparing the results of calculations neglect- 
ing and including 8 respectively. In the same year Brooks (1941) found 
that computed values for the diamagnetic anisotropy of several aromatie 
ring systems are nearly unaffected by the inclusion of S. Recently, 

B. H. Chirgwin and C. A. Coulson (in press), in a more elaborate analysis, 

have proved that mobile bond orders and z-electronic charge densities 
remain completely unchanged, while they concluded that the various 
self- and mutual polarizabilities (introduced by Coulson and Longuet~ 
Higgins 1947a) are only slightly affected. Finally, unpublished 
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calculations by C. A, Coulson indicate, that for a large number of 
polynuclear aromatic hydrocarbons the quotient B/y, y and B being the 
parameters in which the total z-electronic energies are expressed, when 
including and neglecting S respectively, is a fairly constant number. 
The last regularity is of the same character as—and obviously related 
to—the one concerning resonance energies. | 

Only in the case of bond orders and charge densities has a theoretical 
explanation been given. In all other cases the problem is strikingly 
similar: The quantities considered are the sum of certain contributions. 
from the individual molecular orbitals or from pairs of them (e.g. the 
total 7-electronic energy is twice the sum of the energies of the individual 
m.o.’s and the various polarizabilities are built up from contributions of 
pairs of one occupied and one unoccupied m.o.). Now it appears that 
while these separate contributions, being affected by some factor(s): 
containing 8, may and in most cases do change appreciably and irregularly, 
their sum for some reason to be given either changes very little or is. 
merely multiplied by a number, which is approximately the same for all 
molecules. 

In the present investigation the explanation will be proved to lie in 
the fact that S, the value of which is usually taken to be 0-25, only occurs. 
in the summed expressions as $?, which ts obviously a small quantity.. 
First the total z-electronic energies and the resonance energies and then 
the polarizabilities will be discussed (the diamagnetic anisotropy is at 
present not tractable’ by our methods). It will be shown that total 
a-electronic energies are more consistent than resonance energies. Finally,. 
conjugation energies, conjugating powers and their interrelation (Coulson 
and Longuet-Higgins 1948) will also be considered with the same inclusion 
of overlap. It turns out that here again certain conclusions drawn from. 
the simpler theory, neglecting 8, remain valid. 

While the present analysis together with previous investigations referred 
to above, justifies us in neglecting overlap in molecular-orbital calculations. 
on the ground-state of aromatic and olefinic hydrocarbons, this is no 
longer true in the case of heteromolecules. However, it still seems. 
reasonable to assume that if few and small hetero-parameters occur in 
our equations, the inclusion of § will not disturb the consistency of the: 
results essentially. If many and large such parameters are to be introduced 
the value of the entire simple molecular-orbital method as such becomes. 
very doubtful. : 
1. Toe Tora 7-ELEectronic ENERGY. 


If we neglect S, we can write for the energy of the 7’” m.o., ¢,, 


e==h,p. Pi bot ee aera Og 
Consequently the total z-electronic energy is given by 
rh . 
G=a 2 2e,, ahh atone eeeae eee 2-4 (2) 
j=l 


where the summation index, j, runs over the m occupied orbitals. 
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If we include 8, taken to have the same numerical value for all pairs of 
adjacent atoms, we can write for the energy of the r m.o., E,, 


E,.=m,y; eS et Ss) EY, 
where the new energy-parameter, y, is given by 
y=B—Sz«, Meee ee 


‘The relation between m, and k, is known to be (Wheland 1941, Coulson 
and Chirgwin, in press) : 


k 
tes r 5 
Mm, TFS’ (5) 
The total z-electronic energy now becomes 
td . 
é,= pa 2K... . . . . . . . . (6) 


j=1 


‘Our notation implies that «, the Coulomb integral for the carbon centres, 
is taken as the zero level from which all energies are measured. In 
addition it should be remembered that S is a positive quantity, but f and 
y are negative, while & and m are positive for all bonding orbitals. 

From (6), (3) and (1) we derive ‘ 


y Sethe, Bs ti at ; 7 
eked ators ae 


‘This sum can be expressed by the complex integral : 


sey ee CON : 
6-1) las ao | Re a 18 


where the path of integration, I’, is the imaginary axis from — coi to 
+ cov and the infinite semi-circle to the left of the y-axis (fig. 1; compare 
Coulson 1940, Coulson and Longuet-Higgins 1947 a). 4(z) is the secular- 


determinant : 
Ze Pia "Bigs = eae Bin 


Bo —2 Bog. +++ Bon 
A(z)= |... se Ty 


Bra aaa 


Concerning equation (8) we first notice that the pole z= —A/S lies outside 
the contour, Further, if we compare this result with the alternative 
expression given by Coulson and Chirgwin (in press ; See equation (38) of 
their paper), then we observe that in (8) we have got rid of all quantities 
involving Coulson and Chirgwin’s mixed tensor Hy, and instead can deal 
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with the same secular-determinant, (9), as in the simpler theory, neglecting 


8. This is achieved “ at the cost of ” introducing the complicating factor - 
1/(8++Sz) in the integrand. 


If we wish, we may replace (8) by 


‘(z) whe?) 


where » is the number of atoms in the molecule under consideration. 
Replacing (8) by (10) is correct, because 


and it has the advantage that the contribution over the infinite semi-circle- 
vanishes, since for large z, 


Zam Ae) n logge 4 
B+Sz° A(z) B+Sze7S° 2” Sz 
Thus the integral (10) reduces to 
| ay y BY ig B sos | aad PA te: (11) 
ara _oLB+Siy’ Aly) — BF Siy |“ 


Multiplying numerator and denominator of the integrand by (B—Siy) 
yields 


b=t (ee A'(iy) ee aw. (12), 


PLS Ay) RES 


\ 
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43 is an odd function of y, while of course (6?-+S*y?) is even, the 
Ly 


terms 


As 


rain : Aw) and i a 
B+ Sty? A(ry) B-S*y* 
do not contribute to the integral. Hence equation (12) reduces to 
ég= an Es Ay) =F sail lane dy, 
‘ —«@ LB*+S*y*” Aly) B+ Sty” 


7 
which we finally write as 


vou ie 1 neck (gy) 2 
6-35) cee a 


_If we put S equal to zero, we obtain 


ioe) A’ (4 
s=-[" {iu Fay mh ee ME ear 


' which is equivalent to the result of the simpler theory (Coulson and 
Longuet-Higgins 1947 a), provided again that all energies are measured 
from «. 

It will be noticed, that in (13) S only occurs as the small quantity S°. 
Furthermore Coulson has shown (1940), that only the interval 


0<| (y/B)| <4 


need. be considered, in which interval the ‘‘ correction-factor *’ varies from 
i to 3. Finally the integrand in (14) usually reaches a rather steep 
_maximum in the neighbourhood of | y/8 | =2 (Coulson, private communi- 
cation). If this maximum alone determined @, in (13), then we should 
have 


and we can expect the proper result to be only slightly different. In fact, 
as we shall see in the next paragraph, if we equate & and é , as: obtained 
by direct calculations, the result for a large number of molectles is a 
value for B/y of about 0-70—0-75. 


2. RESONANCE ENERGIES AND TOTAL 7-ELECTRONIC 
ENERGIES ; THE QUOTIENT B/y. 

In this paragraph (8/y)¢ stands for the quotient B/y which is found by 
equating & and 6’, while (B/y)y stands for the same quotient when 
obtained by equating the corresponding resonance energies, R and Rg. 
In Table I. values for (B/y) ¢ and (B/y), for a number of molecules and free 
radicals have been collected. Most data for calculating these values have 
been published (see e.g. Wheland 1941, who used B and y in the reverse 
sense); a few new computations were necessary. 
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As already mentioned in the introduction, Wheland (1941) pointed out 
that (8/y), is fairly constant. However, on closer inspection, it cannot be 
denied that rather large differences occur (up to 30 per cent from the mean 
value, 0-50 !), in particular between different types of molecule. On the 


TABLE J- 

Molecule (Bye  (Bive 
Benzene 0-54. 0:73 
Biphenyl 0-51 0-72 
Naphthalene 0-51 0-72 
Anthracene 0-49 0-72 
Phenanthrene 0-50 0-72 
Pyrene 0-49 Oot 
1, 3, 5-triphenylbenzene 0-50 0-72 
Fulvene 0-44 0-73, 
Azulene 0-48 0-72 
Styrene 0-50 0-73 
Stilbene 0-50 0-72 
1, 1-diphenyl-ethylene 0-50 0:72 
1-phenyl-butadiene 0-48 0-73 
2-phenyl-butadiene 0-48 0-73 
Ethylene = 0-80 
Butadiene 0-36 0-75 
Hexatriene 0-39 0-74 
Octatetreene 0-40 0-74 
Decapentzene 0-41 0-73 
Dodecahexzene 0-42 0-73 
Benzyl 0-56 0-72 
a-Naphthylmethyl 0-65 0-72 
B-Naphthylmethyl 0-64 0-72 
m-biphenylmethyl 0-56 0-72 
p-biphenylmethyl 0-58 0-72 
diphenylmethyl 0-53 0-72 
triphenylmethyl 0-51 0-72 
Allyl 0-63 0-75 


other hand a glance at the last column of Table I. shows the values of 
(Bly)g, which were not ‘considered hitherto, to be far more consistent 
(0-73-43 per cent, if for the moment we consider ethylene as a separate 
case). This last fact is in accordance with the theoretical ideas developed 
in paragraph I. 


a a a ee ee eT ee, Ce Oe ne ee ed 
* These hydrocarbons contain odd-membered rings, and are therefore not 
alternant. 
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It can easily be shown, that to a first approximation, a constant value 
of (B/y)g of about 0-73 need not give rise at all to a similar uniformity 
amongst (8/7)g-values. Let us take the case of an alternant hydrocarbon 
with an even number, ”, of carbon atoms. Let us put 


Then according to our assumption above, 


To obtain R and Rg in the usual way, we have to subtract from these 
total 7-electronic energies n/2 times the total z-electronic energy of an 
ethylene molecule (28 and 1-6y respectively). Hence 


R= (K—n)B, 
Rg& (0-73K—0-8n)y. 


Equating R and Rg gives us the result 


(B/vy)pez20-73—0-07 Ey 5 H A = A (15) 
o —1} 
nN 


As K/n is greater than 1, (8/y), can be expected to be smaller than 0-73. 
This indeed is found by direct calculations (Table I.). Furthermore, as 
K/n may show rather large fluctuations, equation (15) cannot be expected 
to yield very similar values of (8/y)g for all types of molecule; e. g. 
comparing benzene with the open chain hydrocarbon with the same n, 
hexatriene, shows us very clearly what happens. For these molecules 
(B/y)@ is 0-73 and 0-74 respectively. For benzene K/n=1-33, for 
hexatriene K/n=1-16. Asaresult we find a large difference in (8/y), (the 
values are 0-54 and 0-39 respectively). In this way, a detailed analysis 
of K/n might give us, for some type of molecule or radical, a more complete 
explanation of the (8/y),z values as obtained by direct calculations, but 
this seems of minor importance to us. 

Obviously the cause of all these difficulties lies in the fact that ethylene, 
which plays such an important part in the calculation of the resonance — 
energies, stands alone in having an exceptionally high value of (B/y)¢, 
namely 0-8. There is no way of getting round this difficulty and the 
important point arising from this discussion is that resonance energies 
calculated with the inclusion of overlap for different types of molecule 
need not be reflected very accurately by the values obtained by means of 
the simpler theory, neglecting S. Total z-electronic energies appear to 
be far more consistent and should be given preference if such is possible 
within the scope of a particular investigation. 
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3. POLARIZABILITIES. 


The definition of the various polarizabilities and several formule to 
calculate them have been given by Coulson and Longuet-Higgins (1947 a) 
within the framework of the simpler theory, neglecting overlap. Coulson 
and Chirgwin (in press) have recently shown how the inclusion of S 
necessitates a change in the formule, while some numerical calculations 
indicate that apart from a factor f/y there is extremely little difference 
between the two sets of figures thus obtainable. 

The short analysis to follow will be given for the specific case of the 
atom-atom “ self-polarizabilities ”’, but it holds equally well for all the 
others. When § is neglected we shall use the customary notation 7, , 
for the self-polarizability of carbon centre 7 in a particular alternant 
hydrocarbon; when § is taken into account, we shall write (7,_ ,)s 
instead. In Table II. a few calculated values have been collected to 
illustrate the similarity between the two cases, for whose existence we must 
find some explanation. 


TaBxeE II. 

Hydrocarbon residue Balint ushty, ea 
Phenyl 0-398 0-359 
2-Butadienyl 0-402 0-377 
Vinyl 0-500 0-469 
1-Butadienyl 0-626 0-601 
1-Hexatrienyl 0-685 0-661 


Without including 8, we have (Coulson and Longuet-Higgins, 1947 a) : 
er 

i ke D a rion 

; j=1 h=m+1 & —E,, 

ie e; stands for the energy of the j’* molecular orbital, as before, and 

; is the coefficient of the r atomic orbital in this m.o. The summation 

lee j and h refer to the occupied and unoccupied molecular orbitals 

respectively. (In (16) the index h is used instead of the more customary 

index k; this is done to prevent confusion with the 4, introduced in 
equation (1).) | 

When including S, equation (16) has to be replaced by (Coulson and 

‘Chirgwin, in press) 


(16) 


4B 9 Cri Crh 
sisleg pee tt (17) 
Mise Y jai homer AGAR(E7—€n) 
~where 
A? =1/(1+k,8) is) 
A= 1/(1+4,8) 
and where k,; and k, are related to «; and e, according to equation (1). 
‘Substituting (18) in (17) eons 


a CF, 
Balke Cee pbk Shh 1 eae 


Y j=1 h=m+1 Te Sh 
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The double sum in the r.h.s. consists of three terms of which the first one 
is evidently B/y.7,,,. The second term, the one involving S to the first 
degree, vanishes owing to the symmetric distribution of occupied and 
unoccupied m.o.’s around the zero-level (Coulson and Longuet-Higgins, 
1947 b). So finally we obtain 


ts. . m n kjk, oe 
Y.(%» g=B Tp ptt? 2 2 » CosCone << eh top 
: 4 j=l h=m+1 k3—k, y 


So the difference between y. (7,,,.)3 and 6 .7,,, will always be small as 
it only involves S?. Furthermore this ‘“ correction-term ”’ will always 
be negative (compare Table IT.) as all terms inthe sum in the r.h.s. of (19) 
are negative. If we wish, we can use equation (19) for the numerical 
calculation of (7,,,)g if 7,,, is known. #.g. for ethylene (one occupied 
and one unoccupied m.o.), with 8 taken to be 0-25 as usual, we obtain. 
immediately 
y(7p,r)s=0°500+4. = ae ; : : ==0-500— = =0-469, 

a result which, of course, is also obtained by the direct application of 
equations (17) and (18). 


A. CONJUGATION ENERGY AND COoNJUGATING POWERS. 


AH. HB AB 


The “* Conjugation energy across a single bond ”’ was defined by Coulson 
and Longuet-Higgins (1948) as follows (see diagram above): ‘ Let us 
connect two alternant hydrocarbons, AH and HB, by eliminating two 
hydrogen atoms at carbon atoms a of AH and b of HB respectively and 
subsequently forming the bond a-b. Then the conjugation energy 
across this bond, R,_,, is defined by (see diagram above) . 


—R,»=¢(AB)—&(AH)—&(HB), . . . . (20) 


where, in accordance with our previous notation, the &’s refer to the total 
7-electronic energies of the species indicated in brackets ”’. 

Coulson and Longuet-Higgins (1948) show that the following relation 
holds in good approximation : 


—R,»/B=constant X B.+/(t1,,4-7,5)) - - » + (21) 


where 7, , and 7, , are the self-polarizabilities of the carbon atoms a 
and 6 in the systems AH and HB respectively. On account of this result 
the self-polarizability is considered as the “ conjugating power” of a 
hydrocarbon residue. 
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If we include the overlap integral S in our calculations, we must 
obviously replace the defining equation (20) by 


—(Ray)s=6s(AB)—6,(AH)—@,(HB). . . . (20a) 


Fig. 2 (a). 


Fig. 2 (6). 
0-6S0 
YV (Ta, alg (Xz .3)s 
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Then it is a natural question to ask whether a relationship still holds 
equivalent to the one embodied in (21), namely 
—(Rap)s/y= constant X y.V/[(Ta,a)s:(7, os] + ~~ (21a) 


If so, then (7,,,)g3 and (7,,,)s act as conjugating powers, just as 7 « 
and 7, , did before. 
22 
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A theoretical treatment is not possible here. It must be remembered 
that equation (21) was never derived rigorously, but even the equivalent 
to the approximate treatment, as given by Coulson and Longuet-Higgins 
(1948) cannot be given here, for the integral expressions for & and the 
(7,,,)s’8 become too involved for this purpose. So we can only check 
equation (21a) by direct numerical calculations. In fig. 2(a) we have 
therefore plotted computed values of —R,,/8 against 8. V/[7a,a 7,01 
and in fig. 2(b) —(Ray)s/y against y . V/[(7a,a)s - (7».)s]. The numerical 
data concerned have been collected in Table ITI. 


TABLE III. 
—Ras/B — B.»/(7a,a-7b,b) —(Rad)s/y —_-y- VL (77a, a)s - (7, b)s 
fe TIES iy 0-383 0-398 0-118 0-359 
eae e ee 
rc 
oP aes o 0-386 0-400 0-121 0-368 
ea hs \ 
C=C 
Cc C 
NY A 
ye Pa “eg 0-388 0-402 hes 0-377 
Yh ie 
C=C C=O 
ESR a 0-424 0-446 0-146 0-410 
4 Dive C )-42 
C 
ne 
roe c=Cc 0-427 0-448 0-147 0-421 
SEE th 
C=C 
ae ieee C=C—C=C 0-460 0-499 0-174 0-464 
ae oy es b 
GHC Get 0-472 0-500 0-175 0-469 - 
a ob 
CH O00... Geet 0-516 0-559 0-212 0-531 
ey fa) 
eto Ue0_ C=C... c=c 0-530 0-586 0-225 0-557 
Coed. 
Gee CeO, 2-3 Cm Om Cac. me nee 0-626 0-263 0-601 
a ob 


A comparison of figs. 2 (a) and 2 (b) shows very clearly that equation (21a) 
holds in as good an approximation as equation (21). This is not a 
complete surprise as we have shown in previous paragraphs that the values 
of & and 7,,, are, in general, consistent with @, and (7, 7)g- These 
quantities completely determine the conjugation energies and the 
conjugating powers. Thus here again the explicit inclusion of the overlap 
integral does not upset the qualitative and semi-quantitative picture as _ 


in Molecular-orbital Calculations for Alternant Hydrocarbons 381 


given by the simpler theory. This result is particularly gratifying since 
some recent investigations show (Coulson, Evans and de Heer, to be 
published) that the conjugating powers in alternant hydrocarbons are 
also of significance in determining certain properties of the z-electronic 
system of substituted hydrocarbons. 


Note added in proof.—Professor J. A. A. Ketelaar and Mr. J. van Dranen 
have drawn the attention of the writer to the following point: It is not 
clear from the manuscript that the complicated analysis of paragraph 1 is 
superfluous if our object is simply to conclude that total z-electronic 
energies are more consistent with respect to the inclusion of overlap than 
resonance energies. They point out that we have, by definition (notation 
as in paragraph 2) : 

Og 4 
Combining the two facts that : 
(i) Rg is small compared with 0-8n and R is small compared with 1-0n- 


(ii) Wheland (1941) has shown by direct calculations that R,/R values 
are fairly constant, 


it is more or less trivial that &;/@ and (f/y) ge are better constants than 


Rg/R and (f/y), respectively. It is also immediately clear from the above 
equation that the (f/y) é will in general be slightly smaller than 0-8 and 


can at most be equal to this value, namely if Rs=R=0 (ethylene). 
In our investigation we have proved theoretically that the consistency 
of &,/€—and thus of (8/y) ge values and their magnitude (--0-8) can be 


expected a priori. We have subsequently shown, that if we take this © 
result as a starting point and then proceed to define resonance energies in 
the usual way, it becomes immediately evident that (8/y)p-values will 
in general not show the same consistency. This second step is entirely 
equivalent to the analysis, given by Ketelaar and van Dranen (see above), 
but in the reverse sense. 
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XXXII. The Estimation of Velocities and Local Mach Numbers in the Flow 
of a Compressible Fluid through Divergent Nozzles, under Adiabatic 
and Isothermal Conditions. 


By R. K. Tempest *. 
[Received January 25, 1950.] 


SUMMARY. : 

In the present paper, a theory is developed giving the approximate 
velocity at any point in an inviscid compressible fluid flowing through a 
divergent nozzle. The method of calculation assumes a radial pattern for 
the stream lines, as suggested by the theoretical and experimental work on 
flows in diffusers carried out by Gourzhienko (1939). The velocity field 
is extended by successive applications of the known solution for radial 
flow from a point source. The formule obtained in the present paper 
are applicable to either subsonic or supersonic flows, they may be applied 
over any length of tube ; and from them the variation of the local speed 
of sound and Mach number may be estimated. 

Numerical examples are given in tables to demonstrate the very large 
increases in local velocity and Mach number in supersonic flows associated 
with comparatively small variation of the expansion ratio of the nozzles. 

A comparison with the results of one-dimensional theory has been made 
in the case where the velocity is approximately uniform across sections 
perpendicular to the axis of the nozzle. It is pointed out that one-dimen- 
sional theory may not be applied to a flow in which a considerable change 
of velocity occurs across a section ; this type of flow is described by the 
theory of this paper. 


§ 1. INTRODUCTION. 
WE consider the flow of a frictionless, compressible fluid in a nozzle with 
axial symmetry. In general discussions of the flow of compressible fluids 
in tubes, linearized equations of motion are used to provide a solution 
valid in some small length of tube. The assumption is made that the 
Mach number and the local speed of sound are constant over this length, or, 
more accurately, variations in these quantities are assumed to be compara- 
tively small and are neglected in the coefficients of terms that are 
themselves small. Numerical evaluation of the velocity at any point is 
cumbersome owing to the complexity of the exact differential equation 
for the potential governing the flow. Under certain circumstances, 
however, these equations may be simplified to give a formula for the 
velocity at any point in very simple form. The theory given below has 
the advantages that it may be applied over any length of tube and that the 


variation of the speed of sound and local Mach number are taken into 
consideration and estimated. 


* Communicated by the Author. 
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We shall first investigate the flow in a cone which is purely radial from 
a source at the vertex and in which the slope of a stream line through any 
point in the field is proportional to the displacement of the point from the 
central stream line. This state of affairs gave Gourzhienko good experi- 
mental agreement in the case of flows which are not completely conical and 
is considered to be a basic characteristic of flows in nozzles. We shall 
assume that in a small element of the nozzle formed by two cross-sections 
at right angles to the axis, this gradation of stream lines exists. Thus the 
flow in such a section will be described by a function which in some 
respects plays the part of a source placed at the vertex of the truncated 
cone formed by the two sections. If the slope of the nozzle at a particular 
section with radius 7p is 69, then we assume that the slope of the stream line 
through the circle of radius r in that section is (see fig. 1) 


(1) 


Variation of direction of stream lines. 


In practice, the use of diffusers has been confined to those of conical or 
almost conical shape, and no experimental evidence on the slope of the 
stream lines in a general divergent flow is available. (Gourzhienko (1939) 
found that in the flow through a cylindrical-conical-cylindrical pipe the 
“assumption of radiality ’ was justified. The flows considered in this 
paper are confined to pipes which have none of the discontinuities in slope 
of the type mentioned above, and we should expect that the assumption 
would be even more justified in these cases. 
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§ 2. NovTaTION. 


We take the axis of the tube as the z-axis ; the coordinate defining 
distance perpendicular to the axis is r, and “corresponding spherical 
coordinates are (R, 6). In irrotational flows, the equation satisfied by the 
potential is 


brit bet b= 3 ($b ret 2b rbobret Bebe.) wt ue ees, (2-1) 


or “: 


t @ 
ss a at serch ae bs 


‘ tes = (#ifant Re PrfoPret pa $s i Poo — ae), - (2.2) 


where @ is the local speed of sound and is a function of velocity if the 
pressure is a function of density only. Equations (2.1) and (2.2) are 
forms of the general equation for the potential in the appropriate coordi- 
nates (Durand 1934). Ifan adiabatic law of the form p=kp” applies, then 


and Bernouilli’s equation assumes the form 


= Il 5 é 
a=" (QP), - - es - + 2.8) 


where q is the velocity at Mach number q/a, and Q is the speed of flow into. 
vacuum, which under adiabatic conditions, is the greatest speed attainable. 
Under isothermal conditions, the speed of sound is constant and equal to: 
a; Bernouilli’s equation is then 


a In p-+4¢=constant. > . 9 0 ae ee 


In this case there is no limit to the speed. 


§ 3. THE FUNDAMENTAL SUPERSONIC FLOW IN A DIVERGENT CONE. 


The solution is exact and the flow is irrotational. The gas state is 
constant on spheres centred at the source, and the velocity increases along 
radial lines from this point. The formula relating velocity to distance R is 
given by the equation of continuity which assumes the form 


pRK*g==codnstant. | .atcg cde ae ene e 
From Bernouilli’s equation (2.3), we immediately obtain 
(Q?—q?)”¥-1gR2=constant. PEs ie ot tine ae 


Let us now suppose that the velocity at a point denoted by R is 
(see fig. 2). 
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It immediately follows from (3.2) that 


Uy [Q?—uP\ Ur-1 R \2 
‘ r a) =(z): Ape ua. = (3.0) 


which gives the outlet velocity in terms of the inlet velocity, the speed of 
flow into vacuum, and the dimensions of the tube. ; 

It should be noted that equation (3.3) holds in both supersonic and 
subsonic cases ; it follows from Bernouilli’s equation that in the super- 
sonic case, Q? is certainly less than (y+ 1/y—1)u?, and for very high speeds 
will be considerably less than this. In subsonic flows Q? is greater than 
(y+1/y—1)ug. 


§ 4. Frow in a Divercent Nozze. 


It is assumed that the equation of the nozzle boundary is r=f(z). 
Generally, the nozzle will be of the form indicated in fig. 3. 


Velocity distribution in a conical flow. 


The initial velocity on a representative stream line is w,; the velocity 
at the point where this stream line intersects a z-constant section is wu. 
It is supposed that the distribution of velocity uw,» across the inlet is known. 
By use of fig. 1, and a knowledge of u as a function of wy, we may find the 
velocity at any point on any streamline. If uw) is constant, the velocity 
will be constant across any section perpendicular to the axis, and we may 
estimate the section at which a particular velocity or a particular Mach 
number is attained. 

We assume that the flow between the two sections z and z+-6z appear to 
diverge from the point on the axis with coordinate z—[f(z)/df/dz], or that 
the slope of the stream lines is proportional to distance from the axis. 
This assumption in effect determines the equations of the stream lines, 
since from (1.1) we may write the differential equation of the family as 


dr AN eaefae die od | 
ae 0 (~) ieee ie hare (In f(z)). Cf. el (4.1) 
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We immediately obtain the stream lines as the family represented by the 
equation 

rapf(Ze <p ls Ce ee (4.2) 


where yz may be any constant between 0 and 1. 


Fig. 3. 


Z+0, 


Velocity distribution in a nozzle. 


Fig. 4. 


We now consider the motion of the fluid contained in the elementary 
volume bounded by the stream surfaces r=pf(z) and r=(u+dp)f(z) 
(see fig. 4). The equation of continuity is now sufficient to determine 
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the velocity wu at any point. The flux across the annulet at z—0 is 
Ser tf (0) } (duf(0)}p9u%) which must equal the final value of this quantity, 
27 {uf (z)} {duf(z)}pu cos «. We therefore. obtain 


(:2 Coe ba 5, 008 ey (Say (Zhe. (4.3) 


In cases of practical interest, the factor {1-++y2(df/dz)? )?}-* may be con- 
sidered to be unity, and the working equation for « assumes the form 
1 


yea ideas Foie 
Q?—u* (9) 

As f(z) > 00, it is immediately seen that w+>Q. It is noted that the 
increase in velocity is determined not by the curvature of the nozzle, but 
only by f(z)/f(0). Thus a short tube with great curvature produces the 
Same increase in velocity as a long tube with small curvature, if the 
ratios of outlet to inlet radii are equal in the two cases. 

It is easy to show that if a divergent flow is being considered 
{Lf(2)/f)] +1}, equation (4.4) has two roots, one lying between u=0 and 
U=U, and one lying between w==u, and w=Q. In the supersonic case, 
wis an increasing function of R and the second solution is the relevant one. 

Exactly the same procedure may be followed to give the subsonic 
velocity of a fluid in a divergent flow. Here the velocity ultimately 
‘decreases to zero and the root of (4.4) lying between u=0 and u=uy is the 
relevant one in this case. 


§ 5. DIVERGENT FLows UNDER ISOTHERMAL CONDITIONS. 


When the equation of state is of the form p=azp, the local speed of 
sound is always da), and Bernouilli’s equation is of the form (2.4). No 
limiting value of the velocity exists as in the adiabatic case. Proceeding 
.as before, we obtain, in the same notation, 


{rat = (B)esr{za e—nd)p Reeraieae? (651) 


In terms of Mach number, 


inay ~ ee? yo (ME Mayes wpe +. «- (8-2) 


§ 6. NUMERICAL RESULTS. 


For purposes of illustration it is now proposed to calculate the velocity w 
at the point at which some representative stream line intersects the 
particular section of the nozzle of radius f(z) in terms of the inlet velocity 
on that stream line at the section of radius f(0). The form of the stream- 
line may be estimated from the considerations of §1. If the streamline 
intersects the inlet section at a distance sf(0) from the axis the equation of 
the streamline is, according to our assumptions, r=<sf(z)/f(0). 
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We suppose that the inlet velocity on the chosen streamline is Up, the: 
velocity on the stream line at the point whose abscissa is 2 is hu, and Q,. 
which is an absolute constant for the whole field of flow, is ku». In super- 
sonic flows, & is always less than (y+ 1/y—1)!, whilst for subsonic flows, 
it is always greater than this value. The inlet local Mach number is. 
obtained from Bernouilli’s equation and is 


Wy h2. \t 
oy 


The numerical work which now arises is the tabulation of a form of- 


equation (4.4) 
1 (WPT \Ur-1 (FIA OS, 
At) {Fat ——€ 4 . . . . . (6.1): 
Fig. 5. 
2:0 


1:5 
f 
Lo 


2 ms 4 

—+> €=f(2)/F (0) 
where « may be interpreted as the expansion ratio of the nozzle with 
reference to section z. If we suppose that the initial velocity u) and 
initial Mach number My are known, then Q=kuy is ealentaeen as. 
{1+[2/(y—1)M@]}!u). We may thus find the value of F(2)/f(0) from 
equation (6.1) which will produce a velocity on that streamline of whup. 
When this velocity is known, the local Mach number is calculated pon 
Bernouilli’s equation. 

The following tables give «—f(z)/f(0) as a function of k and h, assuming 
y=1-4. For purpose of tabulation, if is much simpler to tabulate this. 
quantity than to give a set of solutions of the equation for h when k and e 
are known. From the results of Table I. a graph has been prepared 
giving values of w=Auy as a function of the two variables Q=ku, and 


e=f(2)|f(0) (fig. 5). 


Mach Numbers in the Flow of a Compressible Fluid 389 


The tables are subdivided into four sections involving subsonic and 
supersonic flows under adiabatic and isothermal conditions. As example 
of their use : 

For simplicity let us assume the velocity to be constant across sections 
at right angles to the axis. (If this is not so, we have merely to build up 
the velocity field from considerations of velocities on particular stream- 
lines.) If Q=1-5uo, under adiabatic conditions, the initial Mach number 
(denoted in brackets in Table I.) is 2. In a nozzle of expansion ratio 2-393, 
a velocity 1-3, is attained at the exit, where the Mach number is 3-886. 
Table I a. refers to a similar flow under isothermal conditions. An inlet 
velocity w%) at Mach number 2 is increased to an exit velocity 1-:3u, in a 
nozzle of expansion ratio 1-788. The Mach number at the exit is in this 
case 2°6. 

Thus a comparison is made between isothermal and adiabatic flows at 
the same inlet Mach numbers. 


TABLE I. 
Tabulation of }(z)/f(0)=e (supersonic, adiabatic). Mach numbers are 
denoted in brackets. Results are given to four significant figures. 


1-5 1-75 2-0 2-25 

1-0 1-0 1-0 1-0 1-0 1-0 
(2-981) (2-000) = (1-557) ~~ (1-291) _~—(1-109) 
Ll 1-709 1-200 1-090 1-044 1-019 
(4-141) (2-412) ~—s (1-807) ~— (1-472) (1253) 
1-2 6-112 1-570 1-232 1-113 1-053 
(7-652) (2982) (2-106) _~— (1-677) _— (1-409) 
13 a= 2-393 1 459 1-216 1-107 
(3-886) (2-480) , (1-921) (1 582) 
1-4 = 5-249 1-849 1-369 1-184 
(5-813) (2-982) (2192) _— (1.777) 
15 a 00 2-615 1-602 1-289 
(00) (3-719) (2°535) (2-000) 
1-6 = an 4-613 1-978 1-448 
(5-044) (2-982) (2-261) 
17 es a: 17-40 2-659 1 677 
(9-141) (3-609) —- (2-579) 
18 ne ia =n 4-145 2-030 
(4-617) (2-982) 
1-9 oe a = 9-292 2-622 
(6-990) (3-621) 
2-0 Se a a3 00 3-779 
(00) (4° aa 

= ne co as 6-7 
a (5-812) 
2.2 we. ie = ae 25-47 
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TABLE. LI. 


Tabulation of ¢ (subsonic, adiabatic). Mach numbers given to three 
significant figures. 


3-25 3-5 3-75 4-0 
SNe ie ane oe 
1-0 1-0 1-0 1-0 

(0-723) (0-667) (0-620) _~— (0-577) 

1-068 1-074 1-080 1-085 

(0-568) . (0-525) (0-488) — (0-456) 

1-191 1-204 1-215 1-226 

(0-420) (0-389) (0-362) ~—— (0-339) 

1-423 1-445 1-462 1-477 


(0) (0) (0) (0) 
TABLE La. 


Tabulation of e (supersonic, isothermal, with same initial Mach numbers as in 
; adiabatic case). 


AD es TH 1-0 1:0 1:0 1-0 1-0 
fe (2-0) (1-557) (1-291) (1-109) 
| 1-1 1-520 1-176 1-082 1-041 1-017 

(2-981) (2-2) (1-713) (1-420) (1-220) 

1-2 2-427 1-417 1-192 1-096 1-045 
(3-279) (2-4) (1-868) (1:549) . (1-331) 

1:3 4-066 1-788 1-332 1-169 1-084 
(3-577) (2-6) (2-024) (1-678) (1-442). 
1-4 7-134 2-208 1-513 1-261 1-136 
(3-875) (2:8) (2-180) (1-807) (1-553). 

1:5 13-12 2-849 1-742 1-374 1-199 
(4-173) (3-0) (2-336) (1-937) (1-664). 

1-6 25-30 3-761 2-031 © 1-514 1-277 
(4:472) (3-2) (2-491) (2-066) (1-775). 

17 51-07 5-077 2413 1-686 1-373 
(4770) (3-4) (2-647) (2-195) (1-885). 

1-8 108-1 7-001 2-898 1-895 1-486 
(5-068) (3-6) (2-803) (2-324) (1-996). 

1-9 239-3 9-865 3-530 2-151 1-620 
(5-366) (3-8) (2-958) (2-453) (2-107). 
2-0 554-9 14-20 4-359 2-467 1-780 
(5-664) (4:0) (3-114) (2-582) (2-218). 

21 1347-0 20-88 5-455 2-856 1-972 
(5-962) (4-2) (3-270) (2-711) (2-329). 

2-2 3420-0 31-37 6 916 3-338 2-199 
(6-260) (4:4) (3-425) (2-840) (2-440) 


a i a a ae 
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TABLE Ila. 


Tabulation of ¢ (subsonic, isothermal with same initial Mach numbers). 


h= = 1-0 1-0 1-0 1-0 1-0 1-0 
; (0-791) (0-723) (0-667) (0-620) (0-577) 
| 0-8 | 1-056 1-066 1-073 1-079 1-084. 

(0-633) (0-579) (0-533) (0-495) (0-462). 
0-6 1-168 1-188 1-202 1-214 1-224 
(0-474) (0-434) (0-400) (0-371) (0-346) 
0-4 1-378 1-417 1-440 P4590 el -474 
(0-316) (0-289) (0-267) (0-248) (0-231), 
0-2 1-924 1-972 2-010 2-039 2-064 
(0-158) (0-145) (0-133) (0-124) (0-116), 
0-0 (oe) ee) 0O ee) oe) 
(0) (0) (0) (0) (0) 


§7. PROPAGATION OF CERTAIN VELOCITY PROFILES IN NozzLEs. 


_ A problem of interest to design engineers is the investigation of a non- 

uniform velocity profile in a divergeht nozzle. This problem cannot 
properly be dealt with by one-dimensional theory as the variation of 
velocity across any section may be considerable. The solution may be 
obtained, however, from the Tables I. and II., and fig. 5. A particular 
problem is given below to illustrate the method. 

Let us suppose that at the inlet of the nozzle, radius Ry, we have a 
velocity distribution with magnitudes (1, 0-98, 0-95, 0-87, 0-75)u, at the 
points with radial coordinates (0, 0-25, 0-5, 0-75, 1)Ry, where uy is some 
given velocity. In addition we suppose the flow to be supersonic, with 
a Q-value of 1-5u). The inlet distribution is thus characterized by values 
of & equal to (1:5, 1-53, 1-58, 1-73, 2) at the quarter radii. For the pur- 
poses of the problem, we require the velocity distributions across the 
sections of radii R, and R, where the expansion ratio is 1-5 and 2 respec- 
tively. These may be found as follows. 

At 0-5R, across the inlet, k=1-58. By interpolation from fig. 5, we 
find that if kK=1-58, and e=2, h=1-41. Thus the velocity at 0-5Rs, 
where the expansion ratio is 2, is (0-95) x 1-41), or 1:23u9. By following 
this procedure we may estimate the velocities at the quarter radii of the 
exit sections in terms of uw. For large expansion ratios, we should expect 
the velocity distribution to become uniform and equal to Q=1-5u. The 
following results are obtained, to two places of decimals : 


c=1-0: (1, 0-98, 0:95, 0-87, 0°75)ug at (0, 0-25, 0°5, 0-75, 1)Ry, 
e=1°5: (1.18, 1-18, 1-16, 1-13, 1-10)u at (0, 0-25, 0-5, 0-75, 1)Rs, 
c= 2-0: (1-25, 1-25, 1-24, 1-23, 1-20)uy at (0, 0-25, 0-5, 0-75, 1)Rz. 
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§8. COMPARISON WITH ONE-DIMENSIONAL THEORY. 


We may first recapitulate the results of one-dimensional approximate 
theory relating to flows in nozzles. Whatever the shape of the nozzle, 
-one-dimensional theory assumes the velocity to be constant across sections 
perpendicular to the axis and that this velocity is wholly axial in direction. 
‘This immediately affords the result that 


{7(2)pq=constant;,, a = eee ne 
or, using Bernouilli’s equation, 
{f(z)}*q(Q?—q?)#”-1=constant, .. . . 2 7. (7.2) 


a form which is identical to that of (4.4). Thus, when the velocity is 
-constant across sections of constant z, the theories agree with the proviso 
that the velocity calculated from one-dimensional theory acts at a slope 
O=(r/79)59 Where 5, is the slope of the nozzle at that section. 

In addition, however, the theory outlined in this paper affords values for 
the fluid velocity and Mach number in flows in which the variation of 
velocity across a section is considerable. Such a case cannot be solved by 
the one-dimensional approach. 


§9. COMPARISON WITH EXPERIMENT. 


It is desirable that comparison should be made of the results derived in 
this paper with the results of practical experiment. Application was 
made to the National Physical Laboratory and to the Royal Aircraft 
Establishment, Farnborough, but unfortunately, very little experimental 
evidence is available for comparison. | 

It would also be interesting to compare the results given with those 
obtained from more exact treatment by the method of characteristics and 
it is hoped to carry out such an investigation at a future date. 
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XXXIII. Noise in Current-carrying Ohmic Conductors. 
By Bernard MELTZER", 
[Received December 19, 1949.] 


ABSTRACT, 


The mean-square noise current in a conductor which is in temperature 
equilibrium and carries a current I, is calculated on the assumption that 
the current is due to elementary carriers of charge e. The value obtained 
for the spectral density is {(2kT/R)+[2kT/R)?+ (2eI)?]*}, where R, T, kare 
respectively resistance, absolute temperature and Boltzmann’s constant. 
The interpretation of this result is that, over and above the ordinary 
thermal Nyquist noise, there is present an additional noise, of which the 
spectral density varies in square-law fashion according to (e?R/kT)I2 for 
small I, and in linear fashion according to {(—2kT/R)+2eI} for large I. 
These predictions are compared with Bernamont’s experimental results 
for thin metallic conductors, and the discrepancies are discussed. 


INTRODUCTION. 


In a recent paper (Meltzer 1949) the ordinary thermal noise in a conductor, 
sometimes called Johnson or Nyquist noise, was successfully interpreted 
as pure shot noise. This was done with the aid of an expression, derived 
there, ; 


ela ea... CAT) 


for the average current carried, in either direction, across any cross- 
section of a conductor in thermal equilibrium by the randomly moving 
elementary charge carriers. Here & is Boltzmann’s constant, while T 
and R are the temperature and resistance of the conductor, which has 
its terminals short-circuited ; and it is assumed, for simplicity’s sake, 
that all-the carriers have the same effective charge ¢. (It was shown 
that the result is not changed if the latter assumption is waived.) 
Schottky’s law for pure shot noise, applied to 21), the total average random 
current, immediately gives the correct value for the Nyquist noise. 

The analysis was carried out for the case when there is no source of 
power in the circuit, so that the nett current is zero. Below is presented. 
the generalization of the analysis to the case when the conductor carries: 
a steady current. The result obtained is that in this case, over and above 
the Nyquist noise, there is present additional noise which, for small steady 
currents, varies as the square of the latter, and for large values tends 
asymptotically to a linear law of variation with steady current. 


* Communicated by the Author. 
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This appears to be precisely one of the main features of the extra noise 
observed in thin metallic conductors (Bernamont 1934a, 1934b, 1937), 
and the square-law variation seems to be characteristic also of the noise 
observed in semi-conductors, according to Campbell and Chipman (1949). 
The other chief feature of this additional noise, its variation with frequency, 
is not predicted by the theory in its present state of development, although 
a suggestion is given below of how a fuller theory may account for it. 
A detailed quantitative comparison with Bernamont’s experimental 
results, discussed below, shows that there are some suggestive and 
encouraging agreements with the theory as to order of magnitude, but — 
by no means perfect overall agreement. The discrepancies are discussed. 


ANALYSIS OF FLUCTUATIONS IN A STEADY CURRENT IN A RESISTANCE. 
Fig. 1. 


I 


Battery 


Let the steady current=I. 

If Q is the nett amount of charge carried across any cross-section of 
the conductor in time ¢, the mean square value of Q? is given by a simple 
generalization of Kinstein’s result, quoted in the author’s previous paper : 

Q’= (2kT/R)Yi+ Te. oY Saeivhuse” “eu edaaset oa 

Let us analyse Q, assuming that the elementary carriers all carry a 

charge ¢. We have 


Q=e,+e,+ eee +e,+ eee Ens . . . . . (3) 
where ” is the number of elementary carriers that cross the cross-section 
in time t, and 
the positive or negative sign depending on whether the particular carrier 
concerned was going, say, from left to right, or from right to left. 

' Let p=chance that e;=-+-e, and q=chance that é,=—e, so that 


T=(p—g)(enit)en Ce ee eee (4) 
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Now, from equation (2), 
Q?= Le? 2LEee;. 


Hence at 
Vale + Ize g; 
ids 
fe sn — 1p —g)rete Beene a 3 (3) 
Equating the right-hand sides of equations (2) and (5), gives 
(ne/t)=(2kT/Re)+T?(ne/t)--. 2 2 2... (6) 
Fig. 2. 
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Now (ne/t) is the gross amount of charge carried per unit time across 
the cross-section by all the carriers, those going from left to right and 
those going from right to left. Hence to obtain the noise produced in 
a given frequency band, f to f+-4f, we must apply Schottky’s formula 
(2e Mf times current) to this current. Solving equation (6) gives 

(ne/t)=(kT/Re)+[(AT/Re)?+F]}. 
But we have seen that the amount of current (2kT/Re) accounts for 
Nyquist nois», hence the additional noise is associated with the current 
— (T/Re)-+[(FT/Re)? +P]. 

The application of Schottky’s formula gives the mean-square additional 

noise current 
(i,)P=={--(2kT/R)+[2kT/RP+ (21 4f.  . - . (7) 
22 
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-The graph of spectral density of this noise is shown in fig. 2. 
At first, for small 1, the spectral density follows a square law given by 


(G,/4f=(@R/AT)P2; » . 2. . . ~ (8) 
and then, for large I, it approaches a linear, but not proportional, law 
given by 

”” (G,)/ 4f=— (2kYR)+-2el, . 7 


COMPARISON WITH EXPERIMENTAL RESULTS. 


Bernamont’s first two papers (1934a, 1934b), although historically 
important, may be ignored for our purposes, since—according to his. 
third paper (1937)—the results are vitiated by a failure to take into: 
account appreciable changes of resistance in the platinum wire used. 

In the third paper he dealt with three specimens listed below : 


A: Tungsten (‘‘ Loewe ” resistance) ; 518,000 ohms. 
B: Tungsten (‘‘ Loewe ” resistance) ; 190,000 ohms. 
C: Platinum (‘‘ platine blanc ’”’); 31,200 ohms. 

The results we shall consider are for the total additional mean-square 
noise voltage appearing across the terminals of the specimen over an 
effective frequency band 4f=1200 cycles per second, close to the zero end 
of the frequency scale. The ranges of steady current used and mean- 
square noise voltage observed for the three specimens are as follows : 

-A: 0:3 to 137 microamperes ; 28 to 2-7 10° (microvolts)?. 
7 : 
B: 48 to 490 microamperes ; 9 to 830 (microvolt)?. 


C: 10 to 750 microamperes ; 0-9 to 2700 (microvolt)?. 


If mean-square voltage (v,)’ is plotted against steady current I on linear: 
scales, the curves for all three specimens have the general shape of the 
graph of fig. 2, except that the one for B does not extend into the 
approximately linear region. The asymptotic slopes of the other two 
are given by the following : 

A: 2:36 1014 e.m.u. 
C: 4:75 x 1011 e.m.u. 
Now according to equation (9) of the theory above this slope should 
be given by 
2eR? 4f. 
Hence we may calculate effective values of the elementary charge e : 
A: e¢=0°37X< 10-29 e.m.u. 


C: e=20x 10-29 e.m.u. 
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These are to be compared with the electronic charge of 1-6 10-2%e.m.u. 
If the results be re-plotted on logarithmic scales, one can obtain the 
law of variation over the initial non-linear portions of the curve : 


———9 


A: (v,)?=2-88 X 102°T? e.m.u. 
B: (v,;)?=3-80 x 1047? e.m.u. 
ORY) = 7-94 10? emia: 


According to equation (8), the law should be 
(vpP=(EREAMAT)L ; 
hence effective values of ¢ may again be calculated (assuming T=290 
degrees Kelvin) : 
As e—0-83 x 10? erm.u. 
B: e=0-014x 10-2 e.m.u. 
C 2se=0-30 x 10-*° emu. 
It is seen that for A and C there are discrepancies respectively of factors 
of 2 and 67 in the values of elementary charge effective over the square-law 
and linear regions, respectively, of the observations. On the other hand, 


all the values, except the single one for B, are within an order of 
magnitude of the electronic charge. [Another test of the theory which 


is less encouraging, however, is the following : On the linear scale graphs 


of A and ©, the asymptotes should according to equation (9) cut the 
current axis at [=(kT/Re), so that another estimate of effective e may 
be made from the values of the intercepts. This gives : 


Ae. o==3-4 10 = emiu: 
Geve=-6°8)< 10-2*"elm a. 


It is therefore clear that the analysis given above is inadequate as a 
complete quantitative explanation of Bernamont’s results. Furthermore 
it does not predict any frequency dependence of the noise, which is an 
experimentally determined fact. In fact it seems probable that it is 
just this defect of the theory that is responsible for a large part of the 
discrepancies calculated above. For, some measurements carried out 
by Bernamont on the platinum specimen C showed a drop in average 
spectral density by a factor of 1000 for a 1000-fold increase in frequency ; 
since the effective pass band of the experiments considered above cover 
as much as 1200 cycles per second at the lower end of the frequency scale, 
it is clear that it would have been more appropriate to use results over 
a much smaller frequency band to check the theory. 

A possible explanation of the frequency dependence has been suggested 
by Dr. R. Fiirth in discussion. The effect might be similar to Debye 
and Falkenhagen’s account (1928) of the frequency dependence of the 
conductivity of electrolytes. For the effective charge of an elementary 
particle which crosses a given cross-section may be not merely its own 
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charge e but a value modified by the existence of a charged “ atmosphere ~ 
around it. In electrolytes, such an atmosphere is created because the 
Coulomb forces exerted, say, by a given positive ion produce in its 
neighbourhood an excess of negative ions and a defect of positive ones, 
on the average ; for semi-conductors one might expect something similar, 
electrons taking the place of negative ions and holes the place of positive 
ions; in the case of metals, the main effect is presumably a defect of 
electrons in the neighbourhood of a given electron. If the relaxation 
time associated with the production of such an atmosphere during the 
movement of the particle or hole concerned, is not very much smaller 
than the period of the frequency, one would expect frequency-sensitive 
effects. On the face of it, such effects should produce a lowering of 
effective charge, and it is therefore perhaps significant that of the seven 
values of charge calculated above, all but one are smaller than oie charge 
of an electron. 

The discrepancies described above may also, possibly, be associated 
with experimental error, or the assumption of only one kind of effective 
charge. The theory can be re-stated for a number of different types of 
charge carriers, just as was done for the case I=0 in the author’s previous 
‘paper, but does not then lend itself to comparison with the rather sparse 
experimental data publicly available at present. 


The author wishes to thank Dr. R. Fiirth for stimulating and helpful . 
advice. 
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XXXIV.—A 4-dimensional Generalization of Wilson’s Hypothesis. 


By A. PAaPpaPETROU, 
L.C.I. Research Fellow, University of Manchester *. 


[Received January 31, 1950.] 


SUMMARY. 


Wilson’s hypothesis of electric currents accompanying the rotation 
of electrically neutral matter (in order to explain the origin of the 
geomagnetic field). is generalized in a 4-dimensional language. The 
main consequences of the generalized formulation are discussed. 


SEVERAL attempts at formulating a unified theory of gravitation and 
electromagnetism have been made, but none of them has led to satisfactory — 
results. In spite of this failure, it seems to be almost certain that, if a 
correlation of the two fields really has a physical meaning, it must, first 
of all, explain the origin of the magnetic fields of the earth, the sun and 
other stars in the way recently proposed by Blackett (1949). It therefore 
might be useful to start the discussion from the other end, and try to 
see whether the demand of an explanation of the geomagnetic field could 
give some hints about the nature of the dependence of the electromagnetic 
on the gravitational field. A simple type of correlation of the two fields 
suggested by such a discussion will briefly be described in the following. 
From a phenomenological point of view the geomagnetic field is 
sufficiently described by Wilson’s hypothesis (Wilson 1923). This 
postulates that a rotating material ring will carry the current density 


s=——V(G)pv, Me aoe sani (1) 


p being the mass density, v the peripheral velocity of rotation, and 
G the gravitational constant. Hence it will be sufficient to demand 
from the new postulate, which we are going to formulate, that it must 
reduce to Wilson’s formula (Wilson 1923) when applied to a rotating ring. 
A second demand which we have to impose is that there must be no 
contradiction of the general principles of Maxwell’s theory. In its 
simplest form the last demand is equivalent to the condition that the 
additional “‘ gravitational ” sources of the electromagnetic field, which 
we are going to introduce, must be describable by a charge-current 
4-vector density fulfilling the continuity equation 


Ce ents ea een ea. Fi9) 


* Communicated by Professor L. Rosenfeld. 
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The main difficulty of Wilson’s hypothesis is that it makes use of the 
three-dimensional velocity v. The relativistic point of view demands 
a 4dimensional formulation, and a_ straightforward generalization 
of Wilson’s hypothesis for four dimensions would lead to an additional. 
and extremely strong electric field which does not exist in nature. A 
simple way out of this difficulty is to correlate any given element of mass 
with a corresponding central point, and then to introduce the difference 
of the 4-dimensional velocities of the mass element and the corresponding 
central point into the formula defining the gravitational sources of the 
electromagnetic field. Such a central point is immediately suggested 
in Newton’s theory of gravitation when we consider the gravitational 
field of the earth. The equipotential surfaces const. are almost spheres 
and there is in the interior of the earth one equipotential surface which 
shrinks to a point. We shall call this point the gravitational centre of 
the earth. There is a large number of such centres in space. The moon, . 
the planets, their satellites if they are not too small, the sun, the stars, 
each of them has its own gravitational centre. Any line of force of the 
gravitational field starts from one of these centres and goes to infinity 
without meeting any other centre*. Thus an arbitrary mass element 
is, at a given time, correlated to one gravitational centre; in order to 
find this centre it is sufficient to follow the gravitational line of force 
passing through the actual position of the mass element in the direction 
of the gravitational force acting on the mass element. It is clear that 
this definition of the gravitational centre will be valid in the 4-dimensional 
treatment too, provided that there exists a scalar function corresponding 
to the Newtonian potential of gravitation. This possibility exists even 
in Einstein’s theory of gravitation (general relativity) when interpreted 
as a theory of gravitation in a flat space (Rosen 1940, Papapetrou 1948). 
Besides g;,, which is then only the potential of gravitation, we have 
the metrical tensor ,;, of flat space, and there is the scalar 


b=" (Gix—Nin) 


which differs very little from the Newtonian potential of gravitation. 
(There seems to be no possibility of defining the gravitational centre in 
the usual form of general relativity working with Riemannian space 
only.) 

Having defined the gravitational centre corresponding to a given > 
mass element at a given time, we can now immediately formulate our 
postulate. We call wu! the 4-dimensional velocity of the mass element : 


=| sa} pee 3 
éVO—F)' 1) oe 


Because there are only positive masses. The picture is comparable with 
that of an electrostatic field arising from charges of one sign only. 
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and U' the velocity of the corresponding gravitational centre*. Then 
we put 


AG oleie Ut) ie cise ok tet ton 3 zere{4) 
-and introduce the quantity w fulfilling the equation 
apo Ca tia ae 
Finally we define the gravitational charge-current density by the relation 
a (D=nnp")- eat Ey ee Ab) 
i 


‘One verifies immediately that the vector s; thus defined fulfills the 
continuity equation (2). 

It remains to be proved that in the case of the earth our postulate 
leads to the magnetic field calculated from Wilson’s formula (1). We 
first notice that for all mass elements inside and near the surface of the 
-earth the gravitational centre is a point situated very near to the mass 
centre (or the geometrical centre) of the earth +. Thus if we use the frame 
in which the centre of mass is at rest, we shall have 


Ui=(0, 0, 0, 1). 
We see that, since B=v/c <1, then 
ur V/(l —P?)~1 


-and consequently according to (4) p40 and 
2 I 
pi={p; 0}, p=- V(G)pv, 


v being the velocity of the mass element relative the the mass centre of 
the earth. For a uniform ring around the earth’s axis we shall then have 
. Opt ; 
agi div p=0, 
and consequently from (5) 
w=0. 


** It is easy to see that the position of the gravitational centre will be almost 
simultaneous with the position of the mass element. The gravitational field 
of the planets is nearly static (because of their velocities being small compared 
with c) and consequently the lines of foree—. ¢. the trajectories of the hyper- 
surfaces 6=const—will almost be situated on the hyperplanes ¢=const. 

+ The main displacement of the gravitational centre of the earth relative 
to its mass centre is caused by the sun and it amounts to less than 5 km. For 
the creation of a new gravitational centre near the surface of the earth by 
an accumulation of mass of the earth’s density a body of dimensions not much 
-smaller than those of the earth would be necessary. 
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Tt follows from (6) 
s=— : /(G)pv, sta, 


i.e. in first approximation we have Wilson’s current distribution (1) 
with no charge distribution present. 

From the experimental point of view the case of an isolated body 
near the surface of the earth will be more interesting than that of a 
circular ring, since one might hope to be able to detect and measure 
in the laboratory the magnetic field produced by such a body. We 
therefore calculate the magnetic field produced by an isolated body 
according to postulate (6). Since the velocities will be small compared. 
with c, we omit all time derivatives, and find from (5) 


Aw=— div p. 


This is the same equation as that giving the potential of a dipole 
distribution and its solution is 


| at ae | 


We then see from (6) that the lines of current have the geometrical form. 
of the lines of forces of a dipole distribution whose density is proportional 
to p. From the current density s we have to calculate the magnetic 
potential A by means of the equation 


AA=—4rs, 
and then the magnetic field H= curl A. Therefore 
AH=— 47 curl s=47 curl p. 


The integration of this equation gives 
1 
H= —|- curl p dv, 
which after a partial integration leads to the formula 


1 
H——|—[pr] dv. oS ener 


We see that the magnetic field produced by the mass element contained 
in the volume element dv is identical with that given by Biot—Savart’s 
formula for a current element —p dv. 

Since p is linear in v, it follows that if the body performs a complicated 
motion which can be reduced to a number of simpler ones, the total 
magnetic field produced by the body will be the sum of the magnetic 
fields corresponding to these simpler motions*. The following two. 


$a 

* Strictly speaking the motion of a body in the laboratory will influence the 

position of the gravitational centre of the earth, 7. e. the values of Ui. But 

this will be vanishingly small because of the smallness of the mass of the test 
body relative to the total mass of the earth. 
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examples give characteristic consequences of our postulate (6). A 
projectile moving relative to its terrestrial surrounding with the 
velocity V, while the peripheral velocity of rotation of the earth at the 
place of the experiment is vo, will produce the magnetic field corresponding 
to the velocity v=V+v,. (If, e.g. V=3vy and both V znd » have the 
same direction, the moving projectile will produce a field four times. 
stronger than if it were at rest.) A rotor rotating in the laboratory will 
produce a magnetic field which is the sum of the following two fields : 
a Biot-Savart field corresponding to the translation of the rotor’s mass 
with the velocity vy, and a dipole field corresponding to the rotation 
of the rotor around its axis (the magnetic moment of this dipole being 
given by Blackett’s formula). Since in any experiment we cannot 
obtain velocities much larger than v», we see that the magnetic field 
will always be very small. These consequences are essentially different 
from those following from the generalization of Wilson’s hypothesis. 
suggested by Blackett (1949a): Ascording to the latter there will be no: 
general superposition of the magnetic fields corresponding to different 
motions of a body, and there is the possibility of the production of much 
stronger fields either by suitable motions or by suitable arrangement 
of other masses around the test body. 

More complicated consequences follow from the postulate (6) for certain. 
astronomical cases in which large displacements of the gravitational 
centre occur. For a qualitative discussion it is sufficient to consider 
a planet of radius R and a satellite of radius r at the distance a from the 
planet. Assuming for simplicity that both bodies have the same density, 
we find by an elementary calculation that the satellite will have its own 
gravitational centre only if the distance a is larger than a critical value ao 
given by the relation ; 

. aks 

(a)—1)°= att 

Thus a satellite, moving on an elliptic orbit such that the distance @ 
takes values larger and smaller than a, will have its own gravitational 
centre only on the outer part of its orbit (7. e. the part for which a>d)). 
Hence on the outer part of the orbit the satellite will produce the magnetic: 
field corresponding to its own rotation (around its own gravitational 
centre). On the inner part (a <d)) the gravitational centre of the satellite 
will coincide with that of the planet ; therefore the satellite will also 
produce the magnetic field corresponding to its translational velocity 
round the planet*. 
_ Similar considerations will hold for a double star, though in this case 

further complications will arise if the masses of the two stars are 
approximately equal (and if the two stars are not fully separated). 


Pier 
* This again contrasts with the corresponding consequence of the generalization 
of Wilson’s hypothesis suggested by Blackett (1949), according to which a 
body moving freely in a gravitational field will only produce the magnetic 
field corresponding to its own rotation. 
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Anyway, the general qualitative conclusion from the postulate (6) is that 
large periodic variations of the magnetic field of a double star are to be 
expected, if the orbits are strongly elliptic. In this connection it is 
interesting to recall that Babcock (1948) has observed a star whose 
magnetic field is periodically changing. But very little is known at 
present about the internal structure of this star, so that it is difficult 
to see how this observation could give an argument for or against the 
postulate (6). 

As a final remark we ought to stress that the postulate formulated 
above does not constitute a unified field theory. It might give a 
successful description of the one physical feature of a unified theory, 
namely the production of an electromagnetic field by the gravitational 
field of moving (neutral) bodies. Further we must stress that the 
formulation is not purely differential. It involves the position of the 
gravitational centre of the earth (or generally of the gravitational field 
surrounding the test body), which can only be determined by solving 
the equations of the gravitational field. This is a most unusual feature, 
but we do not think that there are a priori arguments against it. It 
is very difficult to see how one would have to proceed in order to set up 
a complete theory using the postulate (6) ; but some way might be found, 
if only there were some direct experimental evidence in favour of the 
postulate. 
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Further Evidence for the Existence of +-Mesons. 


By J. B. Harpine, 
Imperial College, London *. 


[Received February 15, 1950.] 


We have recently observed an event in Ilford G5 200p “ Nuclear: 
Research *” Emulsions, and this we have to interpret in terms of the 
spontaneous decay of a t-meson, 7. ¢. a particle of mass of approximately 
1000 m (where m is the electronic mass). 

A facsimile drawing of the event is shown. The star centre is 6 mm. 
from the edge of the plate. The angle from the upward vertical to track 
K, which has a length of 740 from the glass to the star centre, is approxi- 

“mately 85° at the emulsion-glass interface. In fig. 2 we have plotted the 
grain density (number of grains per 70) v. range from the star centre of 
tracks a, b, and K, and show for comparison the variation of the mean grain 
density with residual range of thirty protons and ten yz.-mesons. Unfor- 
tunately, track c passes out of the emulsion/air surface after 30 « and has a 
grain density of 29+8-3 (the error being the standard deviation on the 
number of grains counted). 

Track 6, although it showed considerable multiple scattering, remained 
almost in the plane of the emulsion for 1300 before undergoing a large 
angle single scatter. It was therefore possible to determine the product 
momentum x velocity (pf) of the particle which produced this track, by 
measuring this multiple scattering. The mean projected angle of the 
track was measured at intervals of 60, and the difference in angle, «, 
between neighbouring cells thus obtained. The mean of twenty readings 
of « gave ~=1-0°+0-15° and hence we are 99 per cent confident that 
16< pB<46 MeV., the most probable value being 25 MeV. The scattering 
of track a (range from star to glass=420 w) was too small for measure- 
ment, and we are thus confident that the spurious scattering is much less 
than 1/10° per 60. The frequent multiple scattering of track K which 
increased towards the star centre was more characteristic of a typical 
a-meson than of a proton. 

After allowing for the shrinkage of the emulsion during processing (for 
details see Harding, Lattimore, and Perkins 1949) it was found that the 
normal to a plane defined by any two of the tracks, a, b, c, was perpendi- 
cular to the third, within 2°, 7. e. these three tracks were coplanar. Track. 


K did not lie in the same plane. 


* Communicated by Sir George Thomson. 
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Before interpreting these experimental results we must first discuss the 
recording characteristics of the emulsion which was exposed for 85 days 
under three metres of ice. In view of the long exposure, there is a certain 
amount of fading, which, coupled with the low development, makes the 
recording of tracks at minimum ionization uncertain. We do not feel 
justified in comparing grain densities of tracks formed in unrelated events, 
to obtain quantitative information about any particular track. Rather, 
we must analyse this event on its own merits, and merely include data from 
other events in order to show that the fading problem cannot alter our 
conclusions. We have obtained grain density measurements on proton 
and -meson tracks that stop in the emulsion and obtain no significant 
difference in the grain-densities of those coming from the glass and those 


from the surface. Therefore the development (temperature method) of 
the plate is uniform. 


Fig. 2. 
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Both the scattering and grain density of track K decrease with increasing 
distance from the star centre. We thus conclude that particle K is travel- 
ling towards the centre of the star and presumably produces it. The grain 
density measurements indicate that particle K must stop, or at least come 
very nearly to rest before producing the star. The only generally accepted 
‘particles which so behave are 7-mesons which give up their rest mass to 
form stars with excitiation energies of less than 150 MeV. The expected 
grain density-range curve for 7-mesons (as constructed from the curve for 
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y-mesons) is also indicated in fig. 2. If, therefore, particle K is a 7-meson, 
this event must have occurred towards the end of the exposure, for the 
fading must be extremely low as the grain-density of this track is markedly 
above the average grain density that we would expect for 7-mesons. 
Assuming tracks a, b and ¢ to be formed by protons, we find energies of 
~ 300, 110-140, 100-300 MeV. respectively. Thus the total minimum 
energy of 510 MeV. could not be supplied by the rest mass of a 7-meson. 
We cannot exclude the possibility that track a was formed by a fast 
electron, but even then its energy must have been greater than 50 MeV. as 
the scattering was too small to measure. Allowing an improbable 
fluctuation in grain density, track c could possibly be produced by a slow 
energy electron. In this way we could obtain a minimum excitation 
energy of 160 MeV. As yet, no evidence for the production of fast 
electrons following the nuclear capture of a 7-meson has been reported. 
Also it seems improbable that track b could be formed by a proton of 
110 MeV. which would have pB=175 MeV. («=0-15°), the observed value 
is only 25 MeV. This track cannot be produced by an electron of 25 MeV. 
for it should then have minimum grain density, which must be <16, since 
we have observed a track of such a grain density in this event. 

It seems probable that a meson must have produced track b. There is. 
clearly insufficient energy available for a a-meson to produce a 
u.-(216m.) meson (for this is the only lighter one known to exist) and to: 
supply kinetic energy to particle a. In no way then can this event be 
interpreted along orthodox lines of the nuclear capture or of the 
spontaneous decay of a 7-meson. 

The fact that the tracks a, 6 and ¢ are coplanar is very strong evidence 
for the spontaneous decay of a particle at rest. Such a process, namely 
the decay of a r-meson, of mass ~1000m., into a 7~ and probably two other 
7-mesons, has already been published by the Bristol group (Brown et al. 
1949). If we assume K represents the track of a t-meson, the grain density 
measurements would now show that the latent image fading is by no means. 
negligible and hence energies that we infer for particles a, b and c (assumed 
now to be 7-mesons) cannot be very accurate. It would seem, however, 
that their energies should lie in the intervals > 50, 10-20, 5-10 MeV. 
respectively. For pB=25 MeV., a 7-meson has kinetic energy 13 MeV. 
Thus for particle 6, there is good agreement between the energies estimated 
from grain counting and scattering measurements. Resolving moments. 
first perpendicular and then parallel to track a, we find that the energies of 
aand care 50 and 22 MeV. respectively. Allowing 450 MeV. for the rest 
energies of the three 7-mesons, we find the rest energy of the 7-meson to 
be 535 MeV. (=1040 m.). 

If we assume the decay products to be j.-mesons, particle 6 will again. 
have 13 MeV. kinetic energy (for p8=25 MeV.) and the energies of a and 
c are 49 and 21 MeV. respectively. These are still in agreement with the 
energies (now >38, 7-15, 4-30 MeV. for a, b and c) that we infer from 
grain counting. Allowing 330 MeV. for the rest energy of the three 
mesons, the 7-meson would have a rest energy of 413 MeV. (=810 m.). 
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It is seen, therefore, that there is good internal consistency if we interpret 
this event as the decay at rest of a 7-meson of mass about 800m. or 
1000 m., depending on the decay scheme. Since this event is explained 
in enn the same way as was the event of Brown et al. (1949), 
the existence of 7-mesons would seem to be truly confirmed and 
as such should be duly considered in any meson theory of nuclear forces. 
From the theoretical point of view, the range of nuclear forces might well 
correspond to an admixture of 7- and t-mesons. Since the particles are so 
rarely observed (see for example, Franzinetti 1950, in which negative 
evidence for the existence of 7-mesons was obtained), they must be rarely 
produced or have a short lifetime. By now, of the order of 10,000 lightly 
ionizing particles produced in nuclear disintegration must have been 
observed in photographic emulsions and yet no evidence for the decay in 
flight of a t-meson has been obtained. Assuming that 7-mesons are 
produced with kinetic energies about three times their rest mass, we 
conclude that their proper lifetime cannot be much smaller than 10~® secs. 
if they are to account for a large fraction of these fast particles. With 
such a lifetime we would then expect to observe more coming to rest and 
it therefore seems that the hypothesis that they are only rarely produced is 
necessary for agreement with experimental facts. 


I wish to express my appreciation to Professor Sir George Thomson for a 
helpful and interesting discussion on this event. 


Note added in proof.—One further example of a particle coming to rest in 
the emulsion and producing three nearly coplaner lightly ionising particles, 
has just been observed. If the incoming particle were assumed to be a 
a-meson, the fading must be negligible and the three ejected particles, if 
protons, would have a total kinetic energy of >600 MeV. Two of the tracks 
could be formed by fast electrons and only in this case would the energy be 
low enough to be supplied by the rest mass of aw~-meson. The event can 
easily be explained as the spontaneous decay of a 7-meson into three 
lighter mesons with total kinetic energy >70 MeV. or >95 MeV. if the 
decay particles are assumed to be p- or 7-mesons respectively. Slight 
distortion, and the shortness of the tracks, makes a detailed study of this 
‘event impossible. 


Department of Physics, J. B. Harpine. 
Imperial College, London. 
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Intrinsic Block Structure and the Strength of Metals. 


By M. S. Paterson, Ph.D., 


Aeronautical Research Laboratories, Melbourne *. 
[Received February 6, 1950.] 


In his recent paper “On the Strength of Quasi-Isotropic Solids ’’, 
Dr. Fiirth (1949) discusses the idea of an “intrinsic block structure ”’ 
of vital importance in the processes of melting and rupture. He also 
suggests that this block structure may be identical with that assumed 
in Bragg’s theory of the strength of metals and the explanation by 
Wood and co-workers of the X-ray line broadening in cold-worked metals. 
However, from the proposed form of the dependence of the block size 
on temperature (increasing with decreasing temperature), Fiirth’s 
hypothesis would predict a decrease of X-ray line broadening and strength 
with decrease of temperature. This is contrary to experiment, for it 
has been observed that the X-ray line broadening and the stress-strain 
curve for metals are higher at lower temperatures (Wilson and Thomassen 
1934, Hutchison 1948, Paterson and Orowan 1948, Paterson 1950). 
It is therefore doubtful whether the mechanical properties of a metal 
are influenced by an intrinsic block structure such as that proposed by 
Fiirth. 

It should also be made clear that Bragg’s theory deals with the 
resistance of a metal to plastic flow, whereas Fiirth’s deals with its 
rupture. In the case of specimens that ‘‘ neck” during a tensile test, 
the “ ultimate tensile strength ” is a measure of the resistance to plastic 
flow at a certain point on the stress-strain curve (see, for example, 
Orowan 1949). Rupture occurs at a later stage in the test, under stresses 
that cannot be determined from the experiment owing to the complex 
stress state then existing in the necked region of the specimen. 
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XXXVI. Notices of New Books and Periodicals received. 


Theoretical Hydrodynamics. 2nd Edition. By L. M. Minnz-THomson. (Mac- 
millan & Co. [600 pp.] Price £3. 


THE first edition of this work was published in 1938 ; in producing the second 
edition the author has maintained the vector presentation of the subject. 
For the most part, the book is concerned with the theory of incompressible 
inviscid flow ; there is, however, a chapter on viscosity incorporating a section 
on boundary layers, as in the former edition, and in addition in the second 
edition a chapter on compressible flow, subsonic and supersonic. 

The presentation of much of the two-dimensional irrotational flow has been 
modified in the second edition by the use of the “ circle-theorem ”’ which 
though primarily concerned with the effects of introducing a circular cylinder 
in a given flow has, through conformal transformation, repercussions on aerofoil 
problems. The corresponding “ sphere-theorem”’ in three dimensional flow 
is also introduced. The treatment, in the reviewer’s opinion, benefits con- 
siderably by way of systematization by these changes. 

For many problems of inviscid flow vector notation can be used with 
advantage as here, but particularly in such topics as rate-of-strain and stress 
in a viscous fluid its limitations are, in the reviewer’s opinion, only too obvious. 

Le. 


Photoelasticity. By H. T. Jessop and F. C. Harris. [Pp. 184.] (London : 
Cleaver-Hume Press). Price 28s. 


THIS is a timely volume on the principles and methods of photoelasticity 
suitable for the engineer. Since it is written essentially for the experimenter 
the reader only requires an elementary knowledge of mathematics and no 
previous knowledge at all of optical polarization on which adequate chapters. 
are included. There are chapters on reduction and interpretation of results 
and on materials and procedure. The book ends with examples of the use of 
the method and with some appendices of a rather more mathematical nature. 


The Extrapolation, Interpolation, and Smoothing of Stationary Time Series with 
Engineering Applications. By NoRBERT WIENER. [Pp. 163]. (New York: 
Technology Press of M.I.T. and John Wiley & Sons, Inc., 1949). Price $4-00 


‘Tuts important and long-awaited monograph is “an attempt to unite the 
theory and practice of two vitally important fields of work—time series in 
statistics, and communication engineering’. It mainly consists of a basic 
mathematical theory of linear prediction, linear wave filters, and related problems. 
like approximate differentiation in the presence of “noise”, developed by 
Wiener during the war and first published in a restricted report. The present 
publication includes also a useful introductory account of the theoretical 
concepts required, and as appendices two reprinted articles by N. Levinson 
giving a more elementary mathematical treatment, and a numerical table of 
Laguerre functions for use with an orthogonal series solution of the prediction. 
problem. ; 
Engineers may find Wiener’s discussion rather too abstract for immediate 
application, though approximations directed towards the representations of 
ideal solutions by realizable electrical networks are indicated. Mathematicians 
and statisticians will welcome a much more comprehensive account than that 
available in the same author’s Cybernetics, even if the development is in places 
still somewhat elliptical; one problem, at least in many statistical fields of 
application, that of first determining precisely the required correlation functions 
or associated spectra, is not seriously discussed. i Poel Be 


412 Notices of New Books and Periodicals received. 


The Physics of Rubber Elasticity. By L. R. G. TRELOAR. [Pp. vii+255, with 
113 text figures.] (Oxford: Clarendon Press 1949.) Price 21s. 


THIS is one of the series of Monographs on the Physics and Chemistry of Materials 
published by the Oxford University Press. It is designed to give an account 
of recent experimental and theoretical work on the elastic properties of rubber, 
especially the large elastic extensibility of natural and synthetic rubbers and 
the type of molecular structure with which it is associated. ' 

The book opens with an account of the general physical properties of rubber 
and the changes in the thermodynamical functions when a sample of rubber 
is deformed. This thermodynamical analysis provides evidence of the reversible 
nature of the deformation of rubbers, and it gives a firm foundation for the 
statistical developments which follow. These are expounded in the succeeding 
chapters which deal with the elasticity of long chain molecules and of molecular 
networks. Different types of strain are then examined, and it is shown that 
there is good agreement between the experimental results and the statistical 
theory in the region of small strain, but deviations in the region of very 
large strain. This topic is considered further in the closing chapter in which 
an account is given of some of Rivlin’s work in this field. There are also 
chapters on photo-elastic properties, crystallization and mechanical properties, 
and relaxation phenomena . 

The author’s remarks following his equation (2-1) are open to objection. He 
claims that this equation is a definition of internal energy: this claim implies . 
acceptance of the atavistic notion of *“‘ heat ’”’ which is made untenable by the 
experiments on which the first law of thermodynamics is based. The experimental 
basis of the first law is provided by the experiments which Joule described 
in his memoir of 1849 to the Royal Society. These experiments were done 
under adiabatic conditions. They indicate the existence of a function of state 
which is called the internal energy and the changes in which are equal to the 
work done on the system under adiabatic conditions. The author’s equation 
(2:1), which refers to non-adiabatic processes, can then be used to define 
“ quantity of heat’. It is also a pity that the author has not followed the 
recommendations of the International Union of Chemistry (1947) and of Physics 
(1948) in his notation for the thermodynamical functions. 

These blemishes do not, however, detract from the very readable account 
of the properties of natural and synthetic rubbers given in this monograph. 

A. R. MIuuer. 


Twelve Lectures on Theoretical Rheology. By Marxus Reiner, [Pp. 163.] 
(North Holland Publishing Company.) Price 9 Dutch florins. 


PROFESSOR REINER’S well-known “Ten Lectures”? have been revised and 
expanded into the present twelve. The new edition is much more clearly and 
accurately printed than the old. The main change in content is the adoption 
of a rigorous tensor notation for the components of strain, with the use of the 
Kronecker delta to separate the isotropic and deviator components, leading to 
some simplifications and generalization of the formule. The book is an 
excellent brief account of the mathematical basis of rheology. Further 
discussions of the connection between “ solid’? and “ liquid ” viscosities and 
of the dimensional basis of “ consistency variables’ would have been useful, 
and the argument at the top of page 133 seems to me unconvincing. 


F. R. N. NABARRo. 


[The Editors do not hold themselves responsible for the views 
expressed by their correspondents. ] 


